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Abstract. Weakly Interacting Massive Particles are often said to be the best Dark Matter candidates. Studies have 
shown that large Dark Matter-photon or Dark Matter-baryon interactions could be allowed by cosmology. Here we 
address the question of the role of the Dark Matter interactions in more detail to determine at which extent Dark 
Matter has to be necessarily weakly interacting. To this purpose, we compute the collisional damping (and free- 
streaming) scales of generic interacting Dark Matter candidates and investigate the effects on structure formation. 
Our calculations are valid provided the Dark Matter particles have experienced a phase of statistical equilibrium 
at some stage during their evolution. By comparing these damping lengths to the scale of the smallest primordial 
structures known to exist in the Universe, we obtain necessary conditions that any candidate must satisfy. These 
conditions are expressed in terms of the Dark Matter particles' mass and either the total Dark Matter interaction 
rate or the interaction rate of Dark Matter with a specific species. The case of Dark Matter interacting with 
neutrinos or photons is considered in full detail. Our results are valid even for energy dependent cross-sections 
and for any possible initial fluctuations spectrum. We point out the existence of new Dark Matter scenarios and 
exhibit new damping regimes. For example, an interacting candidate may bear a similar damping than that of 
collisionless Warm Dark Matter particles. The main difference is due to the Dark Matter coupling to interacting 
(or even freely-propagating) species. Our approach yields a general classification of Dark Matter candidates which 
extends the definitions of the usual Cold, Warm and Hot Dark Matter scenarios when interactions, weak or strong, 
are considered. 
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1. Introduction. 

The damping of primordial fluct uations has b e en con- 
sidere d more than 30 years ago llSilkl Jl968hM : iMisnerl 
l)l968(l ) when baryonic matter was popular. In the early 
eighties, when the existence of Dark Matter became a 
serious possibility, it was quickly understood that fluc- 
tuations made of Weakly Interacting Massi ve Particles 
will r e nder galaxy for m ation much eas i er llGunn et all 
lll978lk IPeeblesI lll982l) : iBond fc Szafavl lll983ln . If the 
Dark Matter interactions are weak enough and the Dark 
Matter mass heavier than a few keV, one indeed expects 
no collisional damping effects, and no free-streaming on 
scales relevant for galaxy formation. Such particles (e. g. 
neutralinos) are now embedded in the well-known Cold 
Dark Matter (CDM) scenario. As a limiting case, par- 
ticles which marginally comply with the free-streaming 
constraints are usually referred to as Warm Dark Matter 



(WDM ) candidates, the archet ype being l|Davis et alJ 
|l98l|) : iBlumenthal et alJ l)l982|) ) particles with a mass 
of ~ 1 keV whose damping length is about 100 kpc. 

While both the CDM and WDM initial conditions re- 
produce rather well the mass distribution of the observed 
large scale structures, the CDM scenario is in danger t o 
yield an exce s s of small galax i es ( ISchaeffer fc Silkl l)l985(l : 
iMoore et alJ Jl998l) : iMoord lll994^ . The WDM model, 
known to reproduce better the observed faint end of the 
galaxy multi plicity function llSchaeffer fc Silkl jl988l) ). was 
revived {e.g. ISommer-Larsen fc Dolgovl (|XQ0g|^ ^ as a pos- 
sible solution to the discrepancy between observations 
and CDM numerical simulations. However, the collision- 
less WDM scenario does not seem to b e the solution 
to all discrepancies l|Knebe et alJ l|2002() . iBarkana et alJ 
1120011) ). especially if one considers the constraints from the 
late reionizat i on sig n al claimed to be fou nd by WMAP( 
ISpereel et ail (|2003|i ; lYoshida et alJ l)2003|) ). 

Some of the potential problems of the CDM sce- 
nario may possibly be alleviated by invoking astro- 
physic al processes (e.g. IChiu et all ll200llh IStoehr et alJ 
< 2002|) ; iDekel fc Deroil l|2003[) ) or perhaps by redoing 
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data analysis of the dwarf and Low Surface Br i ghtnes s 
rotation curves (e.g. Ivan den Bosch fc Swatersl l|200lj) ; 
Ivan den Bosch et al Also, the D ark halo profile of 

the Milky Way may i Prada e t alj 1120041)) be m ore cuspy 
than previously found l|Binnev fc Evansl l)200l[) ). Should 
the still remaining discrepancy between numerical simu- 
lations and observations be due to the nature of Dark 
Matter or not, a thorough investigation of the Dark Matter 
properties seems worthwhile. 

In particular, one may wonder whether taking into ac- 
count the Dark Matter interactions - generally considered 
to be negligible - would help to solve the Dark Matter 
crisis at small scales. 

This is discussed bv lSpergel fc Steinhardtl (2000) who 
consider strongly self-interacting Dark Matt er to comply 
with t he shape of the galaxy rotation curves l|Moore et alJ 
1999h). Howe v er, a c cording to further investigations 



Yoshida et all ( 2000), Bautz et all ((2002)), this sugges- 



tion, at least in its original form, seems not to be the 
answer. 

One can then address a more general question to deter- 
mine which kind of interaction is allowed for Dark Matter, 
with which particles, and with which strength. 

The original idea that led to the introduction of weakly 
interacting Dark Matter particles (that are now embedded 
in the Cold Dark Matter scenario) was to avoid Silk damp- 
ing as well as prohibitive free-streaming. In the present 
paper, we extend this approach to other Dark Matter in- 
teractions, say interactions of Dark Matter with any other 
species, including itself and seek constraints on the al- 
lowed Dark Matte r particles' parameter space. First in- 
sights were given in lBoehm et alJ I^OOllEooilklChen et alJ 
(2002J). We here undertake an exhaustive and extensive 
analysis by considering also new aspects, considering in 
particular the way the Dark Matter particles achieve their 
relic density. 

We start by relaxing the assumption of collisionless 
Dark Matter to estimate the Dark Matter collisional and 
free-streaming damping effects. Requiring that these two 
damping mechanisms do not prevent the formation of ob- 
served structures allows us to define the range of Dark 
Matter mass and interaction rates that are compatible 
with observations. This is independent of any particle 
physics. We then classify all generic Dark Matter candi- 
dates in a 2-dimensional parameter space (interaction rate 
versus mass) and determine what regions of this parame- 
ter space are allowed by the present observations. 

Our calculations are valid for massive neutrinos, pro- 
vided their mass is within the known experimental and 
observational limits. For the sake of comparison with the 
general case we treat here, we use at various places a 
"reference", the cosmology in which i) the radiation has 
two components (namely photons and massless neutrinos, 
both having the standard temperature) and ii) the matter 
is made of two species, the baryons and the Dark Matter. 

As pointed out in iBoehm et al.l l|200lj) . the collisional 
damping scale of the Dark Matter fluctuations can be ex- 



pressed in terms of two separate contributions, namely the 
self-damping and induced- damping contributions. 

The self-damping is related to the collisional propa- 
gation of the Dark Matter particles. It may be directly 
compared with their free-streaming scale. 

The induced-damping corresponds to the (collisional 
or free-streaming) damping acquired by a species i and 
transmitted to the Dark Matter fluctuations as long as 
the coupling between Dark Matter and this species is large 
enough. The induced-damping effect can be seen as a gen- 
eralization of the Silk damping (of baryons coupled to pho- 
tons), here applied to the case of Dark Matter particles 
coupled to any possible species i. We especially focus our 
attention on the coupling to neutrinos and photons, seen 
to provide a quite large damping and yield rather stringent 
constraints. 

Our results are given in two separate papers. The 
present one (referred to as paper I) is dedicated to the 
calculation of the free-streaming and collisional damping 
scales associated with any interacting Dark Matter parti- 
cles. We then derive the mass and interaction rates of the 
Dark Matter candidates allowed by structure formation. 
This is a rather technical paper. The second one (paper 
II), deals with the physical and astrophysical relevance of 
the candidates whose damping scales, calculated here, are 
not prohibitive. These two papers have been written in a 
self-contained way and can be read independently. 

The present paper is organized as follows: in Section 
121 we recall the expressions of the different transport co- 
efficients for an imperfect fluid made of several species 
(relativistic or not), and take this opportunity to discuss 
issues for which there is still some confusion in the litter- 
ature. In Section we derive the expressions of the self- 
damping and induced-damping lengths. We also point out 
the existence of a new damping effect, that we call mixed- 
damping. The latter is a special case of induced damping. 
It describes the damping acquired by Dark Matter par- 
ticles that would be coupled to a species which is free- 
streaming. In Section^ we derive the constraints on the 
Dark Matter properties from the calculation of the free- 
streaming and self-damping scales. For the specific case of 
Dark Matter coupled to neutrinos or photons, the calcula- 
tion of the induced-damping lengths and the derivation of 
the associated constraints on interaction rates and cross- 
sections are done in Sections El and respectively. 

The definitions of the quantities used in this paper 
are given in Appendix ^ We derive in Appendix [Bl ex- 
pressions for the evolution of the Dark Matter density, 
slightly simplified to be systematically used in our an- 
alytical calculations. The interaction rates appearing in 
the present work are given in the same appendix. The re- 
gions of the Dark Matter parameter space (Dark Matter 
particles' mass and interaction rates) where the various 
damping lengths take different forms are explicitly given 
in Appendix |0 ■ Finally, in Appendix [Dl we display ana- 
lytical expressions for the damping scales and for the as- 
sociated limits on the Dark Matter interaction rates and 
cross-sections. 
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2. Transport coefficients. 

To study the collisional damping of Dark Matter fluctua- 
tions, we first need to determine the transport coefficients 
of an imperfect fluid in statistical 1 equilibrium. The rele- 
vant coefficients are the shear viscosity, the heat conduc- 
tion and the bulk viscosity. They appear in the expression 
of t he energy-momen tum tensor of an imperfect fluid, see 
e.g. IWeinberd l)l97l|b 

Typically, we are lead to evaluate the transport coef- 
ficients of a mixture made of several species i (including 
Dark Matter) . The key which will enable us to obtain con- 
straints resulting from these transport coefficient on any 
kind of Dark Matter will be to write these transport coeffi- 
cients common to all species in the mixture as a sum over 
specific contributions. The latter may be viewed as the 
share of each of the species, contributing to the viscosity 
of the whole fluid. 2 The requirement, in a second stage, 
that each of the contribution in the sum, separately, must 
not be the source of prohibitive damping is then much 
simpler and will yield necessary conditions. 

We proceed now to get the transport coefficients in the 
form we seek. 

2.1. Shear viscosity. 

To derive its expression, we can start from the form given 
bv IChanman fc Cowling! l(l97(lk rj = YliVi valid for a 
composite fluid made of non-relativistic species i in sta- 
tistical equilibrium. The individual contribution rji of each 
species i is given by rji — Pi/Ti where 

r> = £ Ta (f) 

j 

represents the total interaction rate of species i. 

The specific interaction rate of species i with species j 
will be written as 

r i:j = (eru) y . rij , (2) 

where rij is the number density of the species j and (ov)^ 
is the suitably weighted statistical average of the interac- 
tion cross-section between i and j. Under this form, it 
seems natural to expect that this expression can be ex- 
tended to the case of relativistic particles. Therefore one 
can write 

1 = E F> (3) 

i 

1 We use the word thermal equilibrium for particles whose 
temperature and velocities equilibrate at least locally, but 
which are not in chemical equilibrium. Thermodynamical equi- 
librium is meant for particles in both thermal as well as chem- 
ical equilibrium. Statistical equilibrium stands for particles ei- 
ther in thermal or in thermodynamical equilibrium. 

2 It is worth emphasizing that the statistical equilibrium of 
the fluid ensures that each of the species experiences the same 
transport coefficients, given by the sum of all contributions. 



where the sum over i runs over all the species in statistical 
equilibrium with the fluid. The pressure^ may be written 
as Pi = i vf pi, where Vi is the particle's r.m.s. velocity 
in units of c. This yields the shear viscosity coefficient of 
a fluid made of relativistic and non-relativistic species in 
statistical equilibrium: 

One can understand the origin of the expression Q, 
and convince oneself that it is valid also for relativistic 
species, by rederiving it directly from kinetic theory. The 
transport of momentum is induced by a gradient in the 
collective fluid velocity V carried by all species of the 
fluid. The momentum flux induced by this transport of 
momentum may be written ryVV, and corresponds to a 
term rj{dV a / dxb + dV b /dx a — |<5 a {,divV) in the spatial 
part T ab of the energy-momentum tensor. 

Each species i is expected to contribute its own share 
to the total transport of momentum, which we denote as 
rji VV. We may derive the quantity rji by considering the 
momentum flux carried by each of the microscopic par- 
ticles of species i. This momentum flux is v eV where 
v = q/e, q is the particle's velocity, q the particle's mo- 
mentum and e its energy: each particle of species i hence 
carries the collective momentum eV. It is transported 
over a distance of the order of the particle's mean-free- 
path v/Ti, in the direction v, inducing a net transport 
veV(x) — veV(x — v/I\). The latter may be rewritten as 

2 

|^ /Ti VV, taking advantage of the smallness of the mean- 
free-path and the isotropy of the momenta. After summing 
over all momenta of particles of a species i, the flux due to 
a gradient in V is then seen to be ft/I^W, which trans- 
lates into Pl /Ti (dV a /dx b + dV b /dx a - §<5 ah divV) in T ab . 
The coefficient rji associated with species i is thus given 
by Vi — Pi/^i- This is similar to the usual result for the 
shear viscosity of a single species i. The interaction rate 
Ti, however, is the total interaction rate, i.e. it includes 
all possible interactions with other species present in the 
composite fluid. Since the latter arguments do not depend 
on whether species i is relativistic or not, the result holds 
for both relativistic and non-relativistic species. 

The total momentum transported by all species in the 
composite fluid is the sum of the amount of momentum 
transported by each species. Therefore, the shear viscosity 
coefficient of a composite fluid is given by rj — X^ 7 ?^ 
which justifies J3J and (0J. 

2.2. Heat conduction. 

One can also derive the expression of the heat conduction 
coefficient for a composite fluid (including relativistic and 
non-relativistic species) by using kinetic theory. The ex- 
pression of this coefficient for non-relativistic particles has 
not been given in the appropriate form by Chapman and 
Cowling. However, they mentioned that it should be sim- 
ilar to the one obtained for the shear viscosity coefficient. 
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A given species i contributes to the transport of en- 
ergy induced by a gradient of temperature through its 
own conduction coefficient A^. The latter is proportional 
to the particle's mean- free-path Uj/Tj, times a flux factor 
Vidpi/dT. The difference with the previous case (|2.1|l is 
that, here, we consider the energy flux due to a change of 
temperature. Specifically, we get 3 



A,T 



Piv} d In pi 
JY~ d\nT 



(5) 



which holds whether the fluid is relativistic or not. The 
total energy transported by all species is the sum of the 
individual contributions. Yet, the heat conduction coef- 
ficient of a composite fluid is given by A = Y^i^i- This 
finally yields 



AT = 



E 



Pivf d In pi 



3 <91nT 



(6) 



2.3. Bulk viscosity. 



The bulk viscosity is related to departures from equilib- 
rium induced by a non-zero divergence of the velocity field. 
The latter induces a compression or expansion of the fluid 
elements which may result in local conditions differing 
from the global ones. The former have to be compensated 
to recover the general statistical equilibrium conditions. 
This required adjustment provides for another source of 
viscosity. 

Departures from chemical equilibrium also induce bulk 
viscosity contributions. The latter h ave already been ap- 
plied t o other, specific, problems (e.g. lHaensel k, Schaefferl 
( 1992)), and may be included by using similar techniques. 
They allow to constrain the chemical reaction rates, in 
addition to the collision rates we consider in the present 
paper. 

We consider in detail below the departures from ther- 
mal equilibrium, specifically. All the conditions that we 
obtain in this paper on the collision rates thus are neces- 
sary conditions. They must be satisfied separately, keeping 
in mind that there may exist additional constraints on the 
chemical reaction rates. 

It is worth to note, at this stage, that our calculation 
of the bulk viscosity due to departures from thermal equi- 
librium turns out to go much beyond what may be found 
in the existing literature. 

2.3.1. Single fluid. 

A straightfor ward resolution o f the transport equations, 
as outlined in iBernsteinl {HISS}, yields for a single species 



(single) _ dpj/dlliT 



f dpj/dhiT 
\d P i/d\nT 



(7) 



The derivatives with respect of the temperature are to be 
taken at fixed volume. 



The rate Fj is the collision rate for particles of species 
i with themselves. The coefficient u>i is a statistical aver- 
age over momentum-dependent functions, describing how 
the stretching of the momenta by the expansion (or com- 
pression) directly affects pressure and density while the 



ratio 



,/glnT 



dp /dinT describes how the final equilibrium distri- 
bution (i.e. the temperature) is affected. For a single fluid, 
the source of bulk viscosity is therefore the offset between 
the expansion (or compression) which affects the momenta 
producing a non-equilibrium distribution and statistical 
equilibrium which tries to cope with this change by read- 
justing the temperature. The Boltzmann equilibrium dis- 
tribution turns out to be preserved for an ultrarelativistic 
fluid (p — > p/a inducing exactly T — > T / a) as well as for 
a non relativistic fluid of conserved particles (p — > p/a 
inducing exactly T — > T/a 2 ), so the bulk viscosity van- 
ishes in both cases. In intermediate situations, this is not 



so: the deviations of u>i from 



dp, /d In T 
dpi /din T 



however are quite 



small. A numerical estimate shows that (^ sm9le ^ i n g en _ 
eral barely exceeds the 1 % level. This ratio however is 
sizable for the case of zero chemical potential annihilat- 



ing massive particles, for instance Q 



(single 



0.5 at 

T ~ m/20, reaching asymptotically 2/3 for the (unrealis- 
tic: it corresponds to a number density ni oc e~ m l T — > 0) 
case m/T — > oo. 

The result d) for a single flu id is consistent with 
Bernstein's result IBernsteinl ifi988'). However, Bernstein 
considered only massive particles in the limit m/T <C 1, 
that is ultra relativistic particles. The cont r ibutio n 
does not exist in the approach of IWeinherd lll97lh who 
considered either strictly massless or massive non relativis- 
tic, stable (m/T 3> 1) particles. In the latter two cases, 
indeed, the contribution J7J vanishes. As we shall see in 
the next section, there is additional bulk viscosity gener- 
ated in a composite fluid, which is the question addressed 
by We inberg. Despite the comparison made in l)Bernsteinl 
(jl98S^) ^) . these two authors describe different physical phe- 
nomena, with no overlap of their results. 



2.3.2. Composite fluid. 

For a composite fluid, the total bulk viscosity is not simply 
the sum of (|7|) over all species i which are coupled together. 
There is an additional source of inelasticity (compared to 
the single fluid case) which originates from the constraint 
that each of the components i should acquire the general 
temperature of the medium rather than the one it would 
acquire if it was isolated. This implies a rearrangement of 
the equilibrium distribution directly related to exchanges 
among the various species of the fluid. This is in contrast 
with the transport of collective motion or temperature, for 
which each species contributes independently. 

This calculation was first done by IWeiriherel (|l97lh . 
T he latter i s enlightened by an interesting discussion due 
to lZimdah]||)l996|i . 
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With the assumption that the species i are in statisti- 
cal equilibrium, one gets, for a composite fluid : 



( = 5p 
with 



(8) 



5 P -j2 dPl/dlnT 

I 

s p = 2^ — r — 

and 

^d Pl /d\nT ^dpi/dlnT 

w = — f" — w * / L — f" — • 

The weighting by the inverse of the interaction rates {i.e. 
by the collision time of each species i with all species 
within the fluid) is new. It arises from the actual resolution 
of the transport equations in the relaxation time approx- 
imation, the key being to use different collision times for 
each species (this is the very same approximation as in the 
previous cases for the shear viscosity and heat conduction) 
and to account for energy conservation. 

The expression (JSJ may be rearranged as 



i 

with 



(9) 



d Pl /dlnT 



dpi/dhxT 



dpi/dlnT 

d Pl /d\nT f dpi/dWT _ 5_p s 2 
fi \d Pl /d\uT ~ Tp 



(10) 



The first term corresponds to the bulk viscosity of a 
single fluid due to the fact the temperature of species i 
has to adjust to any expansion/compression, which is a 
constraint when particles arc neither fully ultra-relativistic 
nor purely non relativistic. It differs however from J7J) by 
the interaction rate Ti which is now the collision rate of 
particles of species i with all the species in the composite 
fluid (and not only with themselves). Indeed, all collisions 
contribute to the adjustment of the temperature. 

The second term of (|10|) is a contribution specific to a 
composite fluid : 



c 



(composite) dp t /d\nT ( dp t /dh\T _ dp 
Ti \dpi/d\nT Sp 



(11) 



It is due to the fact that the temperature of species i also 
has to adjust to the change of the general temperature 
of the medium under the expansion/compression. Such a 
change represents another constraint when the equations 
of state of the medium (and of species i) are not the same; 
more specifically when g^'yg }"y ^ Indeed, the contri- 
bution |(ri]| may be viewed as the pressure offset due to a 



departure from equilibrium, by A/j = ATj dfi/dT, of the 
Boltzmann distribution /j of species i, in-between two col- 
lisions, with AT, = jjj ( g^g£r - divV. Note that 
this form induces no offset in the total energy-density, as 
required by energy conservation 4 . 

The form 111 111 is in a greement with the expression de- 



rived bv IZimdahll (|199(tI) . However, our result is obtained 
under the simple and unique assumption that all species 
are coupled and in statistical equilibrium. We also take 
into account, in accordance with the style of the present 
paper but to contrast with previous work, that the cou- 
pling of each species i with t he whole medium is, as a rule, 
different. On the other hand, IZimdahll (|l996^ makes much 
more restrictive assumptions, since they are sufficient to 
make his point. He considers species i which are already in 
equilibrium with themselves at different internal tempera- 
tures Ti, under processes he does not show explicitly, and 
uses for his actual calculation a unique interaction rate for 
the va rious spe c ies jus t to equilibrate their temperature. 

In IZimdahll l|199ffl . t he author n oted his result was 
not the one obtained bv IWeinberd l)l97l[) . but did not 
provide the explanation for this. Our form i|ll|l readily 
shows that the contribution of the photons (i = 7) to the 
sumfl specifically is the (sole) term retained bv lWeinberd 
(1971). It is implicit in the latter work that there is a 
second contribution due to the matter component (i = 
to) which however was omitted on grounds the matter 
interaction rate (T m in our notation) is assumed to be very 
large. But due to cancellations in the square of fTTjl. the 
pre-factor does not determine the size of the contribution. 
Indeed, this omitted second contribution is likely to be 
in most cases the dominant one. This is the explanation 
for the difference between these two authors, and will be 
discussed in Section 12 . 3 . 31 using the same composite fluid 
these authors have considered. 

The contribution (|11|) to the bulk viscosity vanishes 
for a single fluid. So, as a rule, in 1)1 ljl. when the coef- 
ficient dpi/d inT j n f rom - Q £ bracket is large, the two 
terms within the bracket nearly cancel each other since 
the composite fluid is dominated by the single compo- 
nent i. Conversely, when the bracket is of order unity, it 
is necessarily for a subdominant species for which there is 
no cancellation. So, C m general never gets very large as 
compared to 77. This is discussed in Section 12.3.31 for the 
generic example of a composite fluid made of matter and 
radiation. In all cases, we find £ < |to if not £ <C |?7. 
As we will see below, in Section |3 the viscosity terms en- 



ter the damping expression with a weig ht C + § »?. For the 
sake of the present study, we simply conclude that the bulk 
viscosity will not induce any new constraint as compared 
to the s hear viscosity , a pro perty which has been already 
used in iBoehm et al.l l)200lh . So, we do not need to con- 
sider the bulk viscosity for the limits on the interaction 
rates we seek. 



4 This can be explicitly checked: Ap = y~) J d 3 pieiA fi — 
^\ ATi/T dpi/d\nT indeed vanishes due to the weighting, 
originally required to this purpose, by the factors 1/IY 
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Our result are far more general than what is needed 
in this paper and more accurate than the ones which may 
be found in the literature. We give indeed the expression 
of the bulk viscosity for a composite fluid as a sum of sev- 
eral well-identified contributions, due to different physical 
phenomena. This sets into perspective the various calcu- 
lations done by the different authors. Also, we take into 
account the difference between the collision rates of vari- 
ous species, a move that is required by our aim to establish 
limits on these rates. We point out several inaccuracies in 
the literature. Assuming an average interaction rate, as 
done in previous work, leaves the latter conclusions un- 
changed. Needless to say, o ur findings are in ag reement 
with the statements made in iBoehrn et alJ (|200lh . 

2.3.3. Mixture of matter and radiation. 

As an illustration of the above results, we discuss here ex- 
plicitly the bulk viscosity of a composite fluid made of ra- 
diation (massless, relativistic, particles with zero chemical 
potential and a unique collision rate T r ) and matter (non 
relativistic particles with a unique rate r m ), assumed to 
form a unique fluid in statistical equilibrium. We discuss 
the various relevant cases, and provide the explanation 
for the different, and contradictory, results found in the 
literature. We hope this will clear out the matters. 

The bulk viscosity for the composite fluid is the sum 

C = Cr+Cm (12) 

and has the same value for the matter as well as the radi- 
ation. Note that - although Cr is the contribution to the 
bulk viscosity from the relativistic species - the bulk vis- 
cosity coefficient of the relativistic species is not Cr but C- 
Similarly, the contribution of the non-relativistic species 
to the bulk viscosity is Cm> but its bulk viscosity of course 
also is C, the two species being part of the same fluid. 

This bulk viscosity may be usefully compared to the 
shear viscosity of the composite fluid : 

V = Vr + Vm , (13) 

related to the pressure p r of the radiation and p rn of the 
matter by r] r = p r /T r and r\ m = p m /T m . 

The case of conserved matter particles. We first consider 
the case of stable matter particles, with a conserved par- 
ticle number. There is then no contribution of the form 
(J7J. From (jnj, we get 

2 



since 



Cr 



dp r /d\nT 



6p 



(14) 



= 3^ 



Vm/& 
Vr + Vm/ 



and 



dp m /d\nT (2 Sp 



(15) 



T^Vm/ 



Vr + W 8 



Sp = \2rj r + -r) m , 



'In 



bp = Arj r 

This yields 
. _ 4 rj r t]„ 



3 rj r + r)n 



(16) 



For r) r > r) m we have Cr ~ V m /^Vr (< Vm) < Vr 
and Cm ~ Vm/G < Vm- One can see in this example that, 
when the bound Q ~ rji is close to be reached, it is for a 
subdominant species. Indeed, we have C °c rj m «C v ~ Vr- 
For the converse, somewhat less realistic assumption rj r <C 
r] m , we would get C ~ Cr oc r/ r <g; r\ ~ r) m . This implies in 
particular C <C ^V- The bulk viscosity is, as a rule, much 
smaller than the shear viscosity. 

Annihilating matter. Here we consider the case of annihi- 
lating matter, when the chemical potential of the massive 
particles also vanishes. More specifically we consider the 
limit of large m/T (having in mind m/T ~ 20) where 
to is the mass of the annihilating particles and T their 
temperature. 

There is in this case a sizeable contribution of the form 
to Cm as discussed in the corresponding section : 



2 

<3»h 



(17) 



Were the matter not coupled to radiation, this would be 
the only contribution to the bulk viscosity, with then C < 
1 77 since C = Cm and 77 = r\ m in this case. 

For matter coupled to radiation, there is also a contri- 
bution of the form ljll|l specific to a composite fluid. We 
2 

get, with b = j^pz, 



Cr 



and 



d Pr /d\nT 



r,. 



Sp 
5p~ 



(18) 



;>)>■ 



bVn 



d Pm /d\nT (T 



Sp 



(19) 



;t>Vr> 



Vr 



Vr + bV7, 



since 



Sp = I2r lr + —rVm , 



TO 
2^2 



TO 



Sp = Ar) r + —rj m . 

This yields for the contribution ifTTjl to lfTU|) : 

4 rj r brj m 



C 



3 n r + bn n 



(20) 
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It is readily seen, for the contributions (jllfl to f, that 
Cr < f^r- One has also £ m < ^br] m , but this is not a useful 
relation since t > 1. A better bound can be obtained by 
noting that, whatever the value of br] m , one has £ m < ^T] r , 
so the damping constraints obtained from £ m do not lead 
to any new constraint once the damping due to rj r is taken 
into account. This is natural since 112011 implies £ < |?y : so, 
when the constraints obtained from 77 are satisfies, there 
are no new ones which appear due to 

The bound ( < |r; still holds when the contribution 
<fTT|) is added to (: the bulk viscosity is not expected to 
bring in any stronger constraint than the shear viscosity. 

2.4. Decoupling and interaction rates. 

In Section |21 we gave the transport coefficients as a sum 
over all species i to which Dark Matter is coupled. The 
coupling condition may be written as Tdm-i > H. More 
generally, a species j is coupled to a species i when 
Tji > H. This condition ensures that the fluctuations of 
the particles j follow those of particles i. More specifi- 
cally this ensures that the fluctuations of particles j will 
be erased at the same scale as those of species i. 

The interaction rates Tdm-i and Ti-dm are in princi- 
ple not equal. There are several reasons for this. First, the 
number of Dark Matter particles and of particles belong- 
ing to species i may be different. Also, the relevant cross- 
section (cru) dm _ i is weighted by the momentum which is 
exchanged during the collision. This transfer may be quite 
different, depending on whether one considers the momen- 
tum transfer to a Dark Matter particle hitting a particle 
of species i or whether it is the converse. 

From the conservation of momentum, the relation be- 
tween those two interaction rates can be written as 



(21) 



with = pi + pi . So, the interaction rates r,j and 
Tji are quite naturally seen to be different. The condi- 
tion Tdm-i > H, for instance, does not necessarily imply 
Ti-dm > H. This leaves also the possibility for Ti-dm < H 
and Tdm-i > H, or more specifically Ti-dm < Ti < 
H < Tdm-i- In this case, when the coupling between Dark 
Matter and species i is efficient enough to modify the prop- 
erties of the Dark Matter fluid but not that of species i, 
one expects the damping of the freely-propagating species 
i to be transmitted to the Dark Matter fluctuations. This, 
of course, occurs only in some range of parameters. 

The relation (|21|l is due to the internal weightings 
within the statistical averages (ov) defined by eq. 10) since 
the interaction probability is weighted by the momentum 
transferred. This is to contrast with what we denote as the 
(av) cross-section which is the usual cross-section which 
counts the number of interactions. The averages (ov)^ 
and (&v)j i in general are not equal : 



h n 3 (Hi 



L n t {av). 



(22) 



i%j irj ■-» \- - I ji -i 

so that Ei (&v)ij ~ Ej (av)j i where Ei is the average 
energy per particle. Whether species i is relativistic or 



not, one can write E ~ 3Tc 2 /w 2 , so for species i and j 
at the same temperature, there is the rule u| (oU)jj ~ 
vf (<yv) In the rest frame of the universe (where the 
blackbody radiation is isotropic), light particles hitting a 
much heavier particle will undergo substantial scattering, 
so each collision is associated to a transfer of momentum, 
and (av) is of the order of (av). On the other hand, heavy 
particles hitting the light ones are barely scattered and 
there is little momentum transfer, so (av) is, as compared 
to (av), down by a factor v 2 /c 2 , much lower than unity. 5 



3. Damping scales associated with Dark Matter 
fluctuations. 

The transport coefficients defined previously are associ- 
ated with dissipative effects. The latter wash out all the 
Dark Matter primordial fluctuations having a size smaller 
than a well-specified scale, that we aim to calculate in this 
section. 

3.1. The collisional damping length. 

Our basic (and sole assumption) is that there exists an 
epoch where the Dark Matter particles are in statistical 
equilibrium; this implies that Td m — J2i Tdm-i > H is, at 
some time, satisfied. Due to the evolution of the Universe, 
one expects Tdm to decrease and eventually to get smaller 
than the Hubble rate at a time tdeefdm.)- 

During all the period where the condition Tdm > H 
is satisfied, Dark Matter is collisional and all the associ- 
ated primordial fluctuations with a size smaller than the 
collisional damping length are expected to be erased. 

This length may b e obtained explic itly in terms of the 
transport coefficients llWeinberd (|197l|)). We write it as 



/ 



c + i 



V 



XT- 



dt. 



(23) 



Here & = V!. ^.. We ha ve adopted the normalization of 
Efst athiou fc Silk ( 1983) so that the mass-scale associated 
with the length l c d is given by M c d = (47r/3) p m l^ d . The 
integral runs over all the period during which Dark Matter 
is collisional. For any reasonable energy-dependence of the 
cross-sections (see discussion in Paper II), the integral l|23|l 
is dominated by late times, so that the value of l c d is 
determined by the Dark Matter decoupling. 

Each transport coefficients 77, AT is given by a sum 
over i. Here i denotes the species to which Dark Matter 

5 The same phase-space arguments obviously hold for the 
interaction rates governing the transport coefficients. This is 
illustrate d by a well-kn own result for non relativistic particles 
(see e. 0. iBalianl il992h 1. The rate (<fu)y ~ (av)\ij = (crv)\ji 
appears in the Lorentz approximation for light particles i hit- 
ting heavy ones j. The rate to be used for brownian motion, 
when one considers heavy particles j moving in the medium 
of the light ones i, is (ov)^ ~ rrn/rrij (av)\ij. Indeed, since 
mivf ~ rrijVj for species at the same temperature, the latter 
cross-section is down by v^/vf, as above. 
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is coupled including Dark Matter itself. A given species 
drops out of the sum at a decoupling time id e c(dm-i)(< 



tdec(dm) 

Dark Matter ceases to be coupled to i. The damping scale 
then reads : 



- 2 E 



(24) 



where i runs over all the species to which Dark Matter is 
coupled (note that for i = dm one has to replace td ec (dm-i) 

by tdec(dm))- 

As discussed in Section |2~31 we can omit the bulk vis- 
cosity coefficient which is either negligible or at most of 
the same efficiency than the shear viscosity coefficient. 
Therefore, it will not provide any new constraints. 

Three different physical mechanisms can be at the 
origin of damping due to interactions, as outlined in 
iBoehm et all l|200l|) : 

1. the damping due to the interaction and displacement 
of the Dark Matter particles within the fluid. This is 
the self-damping. 

2. the damping due to a coupling between Dark Matter 
and a collisional species i which experiences its own 
collisional damping. The Dark Matter then, within 
conditions we will discuss, acquires the damping of this 
species. This is the induced- damping. 

3. the damping due to a coupling between Dark Matter 
and a freely-propagating species i, as discussed in 
Section 12.41 The origin of the damping of species i 
in this case is the usual free-streaming damping. This 
damping is in the present case expected to be trans- 
mitted to Dark Matter . We will call such a damping 
the mixed- damping. It can be interpreted as a part (or 
a particular case) of the induced-damping. 

3.2. Self-damping and induced-damping. 

To get the corresponding damping length, it is conve- 
nient to separate the sum over all species i into a sum 
over the Dark Matter contribution dm and over i ^ dm, 
that is over all species except Dark Matter. Eq.JSHJ) can 
then be split into two parts. One which only depends on 
the Dark Matter properties (corresponding to the self- 
damping length) while the other one depends on the 
characteristics of species i ^ dm (corresponding to the 
induced-damping length). This yields 



led — Isd + lid 
i^dm 



(25) 



noting that the damping contributions combine quadrat- 
ically. Using the explicit expressions of r)i and A^T given 
by eq.(@J and eq.@, one finds 



3 Jo f>a z ldm, t 



(26) 



and for each species i ^ dm : 



). This time corresponds to the epoch at which ft 



2^2 
3 



i>a 2 Y, L 



(1 + Qi 



dt 

T 



The factor 
P 







2 d In pi 



(27) 



(28) 



obtained by using © and lt2*3l is associated with the 
heat conduction coefficient. Here, it has been calculated 
neglecting the p ressure of no n -relat ivistic particles. As 
already seen in l|Boehm et alJ I^OOllu . this yields A,- oc 
dTnT~ = dWr^ ~ ® f° r a non - reia tivistic species j. The 
heat conduction coefficient is then non-zero only for 6 a 

2 

relativistic species i (&i cx ~2 m ~) an d is important (i.e. 

Qi larger than unity) only in the matter-dominated era 
where we have for any relativistic species i the relation 
Qi ~ esn > i and ~ < 1. 

Pr Pm Pr — 

It turns out to be quite useful - and sufficient within 
the accuracy of our calculations - to approximate the in- 
tegrals (|26|l and H27|) . written as J by the maximum 
the factor K reaches within the integration interval: 



J Kj ~ M&x[K] 



(29) 



It is then straightforward to see that these integrals are 
dominated by the late times. Indeed, the decoupling cor- 
responds to an epoch where the system passes from a cou- 
pled to an uncoupled regime while the Universe is expand- 
ing. This implies that if is a growing function of time, as 
can be readily inferred from a short calculation counting 
the powers of t. 

The approximation (|29|) . which aims to compare pro- 
cesses that differ by orders of magnitude, is good within 
factors of order unity. It considerably simplifies, at a very 
low cost, the overall classification of all kinds of Dark 
Matter we undertake in this paper. 

As a result, the self-damping scale may be written as 



2 2"7T 2 PdmVd m t 
l sd ~ — j_ 2 -n V 1 + ® dm > 



ia 2 T dr , 



dec(dm) : 



while the induced-damping due to species i reads 
, 2 2tt 2 pivft 



3 i>a?Y, 



(i + eoi dec(dm-i) 



(30) 



(31) 



It is convenient to rewrite these approximate expressions 
by introducing explicitly the Hubble rate (written as H = 
a/a = aft, with a = 1/2 or 2/3 depending on whether 
the Universe is in its radiation or matter-dominated era). 
This yields 



hd ~ TTdri 



H 

r flirt 



Vdmt, 



\dec(dm) i 



(32) 



6 The complete expression including the pressure o f the mat- 
ter co mponent can be derived along the lines of IWeinberel 
but it is not needed in the present paper. 
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and 

hd ~ 



H 



\dec(dm—i) 



(33) 



where we have kept the factor H/T (although it is unity 
at the Dark Matter decoupling) in order to extend it to 
the specific situation discussed in section EL 2. 31 Here, we 
have introduced the dimcnsionlcss ratios 



Tdr, 



and 



3a 



(1 + 9dm) 



2ft 
Sad 1 



(34) 



(35) 



From the relation Pi< $> and the discussion in section 1X11 
of the possible values that the factors 6<j m , 6^ may take, 
one can readily see that r^m and ri are smaller or equal 
to unity. 

Note that our formulas are valid for energy-dependent 
cross-sections, that is for temperature-dependent averages 
(av). 

3.3. Free-streaming. 

After its decoupling at time td ec (i), species i is freely- 
propagating. Its primordial fluctuations experience the so- 
called free-streaming damping. The damping acquired at 
a given time t is proportional to the average distance trav- 
eled by a particle of the species, that is in comoving units : 



MO 



V; dt 



(36) 



It has been argued, and checked numerically, that an in- 
ter esting approximation t o this damping length is given 
by llBond & Szalavl Jl983l) 1 : 



MO 



7r Max 



Vit 



(37) 



where Max denotes the maximum value of the free- 
streaming scale within the integration interval \td ec (i)> t] 
of ea. (|36|l . This is a matter of convention about what one 
calls free-streaming length. In the regime where we write 
our constraints, the two approximations are not dramat- 
ically different, and one might as well use the simpler 
one H37fl . With the norm alizatio n 1371) . th e free -streaming 
length lf s is associated l|Bond fc Szalavl 1 19831) ) with the 
damping of all the fluctuations of mass scale 



M 



(47T/3) Pr 



I 3 
Jfa ■ 



(38) 



The quite simple approximation l|37|l is well-adapted for 
the comparison with the collisional damping length which 
is also approximately evaluated by means of a similar sim- 
plification (J2UJ. 

As a rule, close to the parameter values (mass and 
interaction rate) adapted to our constraints, (|3"rJjl yields 



length scales measured in the same dimensioned factors 
as the ones obtained from Ij37(l . but with numerical co- 
efficients larger by a few units due to logarithmic con- 
tributions (but then it is not clear whether the mass- 
scale/length-scale relation is still (|38|) or the one where 
lf s is to be divided by 2). For the sake of comparison, we 
give in Appendix ID . II the expression of the free-streaming 
lengths calculated by means of the approximation (J3EJ. 
The latter might be more suited than (|37|l far from the 
parameter values we handle in this work. 

3.4. Mixed-damping. 

A freely-propagating species may transfer its own damping 
to the Dark Matter fluctuations if the condition Ydm-i > 
H > Ti > Ti-dm is satisfied. This in particular requires 
that 



> 



dm 1 



due to the symmetries of Ti-j under the 
permutation of i and j, as can be inferred from Eq. (|21(l . 
One can then decompose the induced-damping scale into 
two parts. 



/ 



2tt 2 
3 

2tt2 
3 



(1 + 6*) dt 



ft v i 
l/>a 2 H 



(I + 6i) dt 



(39) 



The first term, valid as long as Ti > H, up to the 
time at which i ceases to be collisional, is the (collisional) 
damping acquired by i which is communicated to the Dark 
Matter . The second one is an estimate of the damping 
acquired by i during the period where Dark Matter re- 
mains coupled to this species, now freely-propagating since 
Ti < H. It is obtained by simply replacing 1^ by H to get 
an expression where the distance covered by a particle of 
species i while it was collisional is replaced by the distance 
it covers when being free-streaming. The mixed-damping 
length may finally be written as: 



I 



2^2 
3 



tiecldm-i) p. y 2 



i>a 2 E 



(I + 0i) dt 



(40) 



With the same approximation as in Eqs. ll-i-il) or H37fl . we 
get: 



Vit 

\dec(dr, 



o • 



(41) 



Omitting the factor r^, one easily recognizes the free- 
streaming scale of the species i indicating that the damp- 
ing acquired by i is indeed transferred to the Dark Matter. 
In fact, because we are considering a composite fluid, it is 
not the free-streaming of species i which is communicated 
to the Dark Matter fluctuations but the free-streaming 
weighted by the factor r*j. This length finally contributes 
quadratically to the total damping scale l c d- 

One may however consider that (|41|l is a too naive 
estimate. It is then worth to note that a lower bound to 
the damping length may be obtained by only considering 
the damping acquired by the Dark Matter fluctuations 
during the collisional regime of species i. This contribution 
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corresponds to the first term in l|39|) and represents a lower 
bound to the true damping length since 



& > 



2tt 2 



if) a 2 r 4 



(l + 9i) dt 



(42) 



This integral may be written, with the same approxima- 
tion as in 13. II : 



lid > 7T Ti \decU) 



Such a result solely relies on the standard physics involved 
in the collisional regime. It differs from the mixed-damping 
estimate, ea. (14H . just by the time at which the factors 
are evaluated. The collisional bound is estimated at t dec ^ 
at which the species i enters the free-streaming regime 



while the mixed-damping is estimated at t, 



dec(dm-i) 



the 



true epoch at which the Dark Matter is no longer influ- 
enced by i. This collisional lower bound may be much 
smaller than the actual contribution, that we expect to 
be closer to the mixed-damping estimate. The latter will 
provide substantially stronger constraints. So, for our es- 
timates we will use only this length, although it relies on 
less conventional physics. 

3.5. Dark Matter damping due to its own 
free-streaming. 

After its total decoupling, Dark Matter freely propa- 
gates erasing all fluctuations having a size smaller than 
its own free-streaming length. From |37|). this is at time 



t > t. 



dec(dm) 



lf s {dm) ^TrMax 



(44) 



where the maximum is to be taken within the time interval 
[tdec(dm),t\ where Dark Matter is freely-propagating. 

3.6. Comparison of the free-streaming scales with the 
self-damping and induced-damping scales. 

There is a striking similarity between the expressions of 
the free-streaming scale Q370 and the self-damping scale 
<|32|) . This is simple to understand. Since the collisional 
damping scale is based on diffusion processes, the associ- 
ated damping length is proportional to the distance trav- 
eled by collisional particles during their random walk. 
In a non expanding Universe, this distance is usually 
much smaller than the distance traveled during the free- 
propagation. However in the present case, due to the ex- 
pansion of the Universe, it gets increasingly large up to 
the decoupling time. The collisional damping length is 
hence dominated by the late times, where the random walk 
marginally corresponds to the free propagation. 

The self-damping length loL'l) is however still somewhat 
smaller than the Dark Matter free-streaming length (|37|l 
since the transport of heat or momentum is done with an 



efficiency translated in terms of the factor r^m (smaller or 
equal to unity). We therefore have the relation 



Isd < I 



fs(dn 



(45) 



The factor being also smaller or equal to unity, we in- 
fer that the induced collisional damping due to species i 
satisfies 



(43) l *d < lfs(i) 



(46) 



where lf s {i) is the free-streaming length (|37|l of the species 
i. Therefore, the induced-damping scale of Dark Matter 
fluctuations is not bounded by the distance traveled by 
the Dark Matter particles but rather by the distance cov- 
ered by the most rapid particles to which Dark Matter is 
coupled. This is very important for massive Dark Matter 
particles in statistical equilibrium with relativistic species 
since, in this case, the Dark Matter induced-damping scale 
is bounded by ct and not by VdmX- In this case, although 
the thermal velocity of such Dark Matter particles is rel- 
atively small, their coupling to relativistic species allows 
them to undergo substantial damping. 

4. Constraining the Dark Matter properties from 
free-streaming and self-damping. 

The self-damping and free-streaming lengths happen to 
depend only on two physical Dark Matter parameters: the 
Dark Matter particles' mass rridm and the Dark Matter in- 
teraction rate with the medium Tdm- The dependence of 
the free-streaming damping length on the Dark Matter 
particles' mass is well-known, its dependence on the inter- 
action rate is less obvious and is in some sense the new 
feature examined in this section. 

The constraints we get in this section are obtained 
(there is one exception discussed later) at the time td ec (dm) 
where Dark Matter decouples to any species i, including 
itself. The rate Tdm depends among others on the number 
density of the particles in the medium. Due to the ex- 
pansion, this density, and whence the rate, decreases with 
increasing decoupling time, that is with increasing interac- 
tion strength. It turns out to be convenient to work with 
an interaction rate which, taken at this epoch of refer- 
ence, is an increasing function of the corresponding cross- 
section. Therefore we will use (see also Section lA.2.6|l the 
" reduced" interaction rate Tdm = Td m a , where the most 
obvious time-dependence due to the number density evo- 
lution is removed. 

4.1. Classification of the Dark Matter self-damping 
processes. 

We subdivide the plane (md m ,Td m ) into regions within 
which the free-streaming and self-damping lengths have 
different expressions. The relevant time-scales which de- 
termine the damping lengths of the Dark Matter primor- 
dial fluctuations are 
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— the Dark Matter decoupling time t dec ^ dm j (at which 
Dark Matter has decoupled from all species) 

— the epoch at which Dark Matter becomes non- 
relativistic t nr 

— the usual epoch of equality t eq ^ 1+u y 

These times may be translated in terms of scale-factors 

a dec{dm)i a nr and 0, e q(^+u) = Pj+v (To)/p m (T ). Their 

relative ordering defines six regions in the Dark Matter 
[mass/interaction rate] parameter space: 



Region I 


a dec(dm) < a nr < 0>eq(j+u) 


Region II 


CL nr < Cldec(dm) ^ ®eq(~f+is) 


Region III 


d nr <C deq(~f+is) *^ Q'dec(dm) 


Region IV 


a dec(dm) ^ a eq(-)+v) ^ a nr 


Region V 


a eq(l+v) ^ a dec(dm) a nr 


Region VI 


a eq(Sf+v) a nr < a dec(dm) 



This classification provides a natural way to separate 
heavy from light, late from early decoupling Dark Matter 
particles. This also allows one to distinguish particles 
which are able to annihilate and those which are not. 
Region I includes in particular the standard Hot Dark 
Matter (HDM) scenario (like for instance neutrino Dark 
Matter ). Region II, includes, among others, annihilat- 
ing particles (e. g. Weakly Interacting Massive Particles). 
Region III is devoted to strongly interacting particles. The 
last three regions, for which a eq r 1+JJ \ < a nr , correspond 
to the less likely very light Dark Matter particles, with 
masses below a few eV. 

These Regions are shown in Fig. ^ Explicit equations 
for their borderlines are given in Appendix lC.il We now 
give the expressions of the damping lengths in these re- 
gions, and the associated constraints, whose values also are 
displayed in Fig.^ Analytical expressions of the damping 
lengths may be found in Appendix lD.il 



4.2. Self-damping and free-streaming lengths. 

4.2.1. Region I : a dec ^ dm) < a nr < a eqiri+u) 

Since CLdecVdm) corresponds to the total decoupling of Dark 
Matter , annihilation processes by definition are frozen 
afterwards. The Dark Matter comoving number density 
then remains constant : n,dm(T nr ) a\ r — nd m {To). The 
former can then be translated in terms of energy den- 
sity p dm — m dm n dm , which, at the non-relativistic tran- 
sition (a = a nr ), is required to be equal to the energy- 
density 7 in the relativistic regime : p dm = ^egd m dBTd m = 
\g*draa,BT^ ■ This turns out to be a very stringent require- 
ment. To be satisfied, it requires a quite small effective 
g*dm to ensure that the number density of Dark Matter 



The factors entering the energy-densities are defined pre- 
cisely in Appendix 1X1 . 
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Fig. 1. Bounds in the [Dark Matter particles' mass, 
Dark Matter interaction rate] parameter space obtained 
from^self-damping and free-streaming. The interaction 
rate r<j ec (dm) is ^dec(dm)0? taken at the Dark Matter de- 
coupling. The labels I to VI correspond to different sce- 
narios. The hatches indicate which part of the parameter 
space is forbidden. It corresponds to the regions where the 
free-streaming and self-damping scales are above lOOfcpc 
(i.e. ~ 10 8 M©) scale. This is seen to provide a constraint 
involving both the interaction rate and the mass. The in- 
dications in this figure are schematic. Some factors in gen- 
eral of order unity are omitted for simplicity. The reader 
interested by these constraints should use the expressions 
in the text, where all factors are given explicitly. 

particles at the present epoch matches the observed relic 
abundance. Specifically 



Pj ) 



3.36- 



« 



eq(-y+i>) 



(47) 



This condition may be rewritten in terms of /tj m , the 
factor which determines the ratio of the photon to Dark 



Matter temperature, as Kd m (T nr j 
is 



r 9dn 



3.36 a nr K dm (T nr ) 



that 



9dr, 



0.3 



(^y. (48) 



V \MeV > 



In this scenario, the Dark Matter number density is 
to be set at a quite remote epoch, and all changes in 
the particle number are frozen out, much before the non- 
relativistic transition. We will thus call this way of match- 
ing the present relic density the Ultra-Relativistic Freeze- 
Out (URFO) scenario. 
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The requirement (|48|l may be achieved through vari- 
ous ways. One one hand, Kd m may be large (quite larger 
than unity) if the decoupling between the Dark Matter 
(which may remain interacting) and the thermal bath oc- 
curs at a remote enough time, generically much before its 
total decoupling. On the other hand, one can have a small 
e if the chemical reactions among Dark Matter particles 
which do not conserve the Dark Matter number density 
freeze out at early enough an epoch, the contact with the 
thermal bath being maintained. The particle number con- 
servation from this epoch till nowadays then ensures that 
the Dark Matter density has the appropriate value at the 
non-relativistic transition and at the present epoch. 

The expressions of the self-damping and free-streaming 
scales in this region read respectively : 



sd 



ct 



r dm \dec(dm) 



and 



ct 



fs 



(49) 



(50) 



Numerically, one finds the results given in Tabled where 



r 



dec(dm) 



is the reduced Dark Matter interaction rate at 



td. pr(d.m.\ ■ The free-str e aming scale corresp onds to the result 
of lBond et all lll980h : iDavis et alJ lll98lft . 



Table 1. 



1 dec(dm) 



ll J = 0.35 kpc 

c(dm)) Q * (^dec(dm)) ^q_29 s — 1 

if, = 0.51 kpcg'-J{T n r) (™#^) _1 



The comparison of l|49|l with (|50|l shows that the self- 
damping contribution is down by a factor 



ct I 

a \dec(dm) ^dec^dm) 



< 1. 



(51) 



The self-damping is therefore generically smaller than 
the free-streaming in Region I and will not be 
the most constraining effect. The coefficient r^m ~ 

)/^(7dec(dm))] ^ j is also below unity as long 
as Dark Matter is coupled to radiation. Indeed (as may 
be inferred more readily using the ratio of the energy- 
densities given in Sect. IB.2|l pdmj V ~ Pdm/Pr fs a- nr ja eq 
implies rdm fs 1- On the other hand, if the decoupling from 
radiation occurs before the Dark Matter total decoupling, 
the Dark Matter fluid does not contain any radiation com- 
ponent and one gets: $ ~ p rn and tdm ~ 1- 
From (|50[1 . we get the usual condition 



(T nr ) >5.1 keVg'^(T nr ) 



( Is 



1 lOOfcpc J 



(52) 



Implementing the estimated value l|48(l of K& m , yields the 
condition 



mdm > 0.9 keV 



( 9dm\ 4 ,- 



O h Z \ 4 / / \^ 
iii m /t 7 n \ / ^struct \ 



(53) 



at the orig in of the so-called " IkeV" mass limit 

l|Davis et al.l <|l98lh ) 



0.3 / \l00kpcj 

t 
))• 



4.2.2. Region II : a nr < a dec(dm) < a eq ^ +u y 

The Dark Matter (total) decoupling in Region II occurs 
after the non-relativistic transition. There are in this case 
two ways of satisfying the relic density requirement : 

— The first one, as in Region I, is to consider the URFO 
scenario where the constraint l)48|) has to be imposed. 

— The second one is to consider particles which remain in 
chemical equilibrium up to their non-relativistic tran- 
sition (actually slightly after, as we shall see below). 
The allowance for this second possibility opens a new 
window. Indeed, one may have a large energy density 
before the non-relativistic transition. Then the Dark 
Matter particles, still in chemical equilibrium, annihi- 
late and the number density exponentially decreases. 
Once it reaches the appropriate value, one assumes 
that the reactions freeze-out. We will call this way of 
matching the present relic density the Non-Relativistic 
Freeze-Out (NRFO) scenario. 

In the NRFO scenario, to match the observed relic 
density, a specific value of the freeze-out temperature, 
and in turn of t he a nnihilation cross -section (see e. g. 
iLee fc Weinberd l|l977(k iGriestl l)l98S|n is required. We 
only shortly discuss its derivation, which is well-known. 
The conservation of the Dark Matter number density in 
the NRFO scenario is relevant only after the chemical de- 
coupling. By equating the comoving Dark Matter number 
density during its exponential decrease (just before freeze- 
out) to the comoving density just after freeze-out, with 



rpfo TT^dm 

dm ~ ~~~ ' 

one sees that x f a obeys the equation 

■ / tic \ 

(2tt 



x 2 fo e ^ 



^ <! in 



T J g*dm/ K 3 dm ( T fo) 

yielding the very simple rule 



X fo 



14 + In 



(Tfo) 
lMeV 



(54) 



(55) 



(56) 



where we have neglected unimportant logarithmic factors. 
For ~ 100 GeV particles, Xf a is of the order of ~ 20: the 
annihilation reactions freeze out somewhat after the non- 



relativistic transition, at T 



Jm = ^ao*- This implies a well- 
defined value for the cross-sections to maintain the chem- 
ical equilibrium up to i = Xf and not after. It yields 
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the standard constraint on the Dark Matter parameters, 
which in this case must be satisfied instead of JJH}. Larger 
values of (cv) ann are prohibited 8 , as well as lower values, 
unless they are small enough so as to comply with the 
URFO scenario and therefore satisfy condition l|48[) . In 
the special case where the chemical equilibrium is reached 
via direct annihilation, for instance, the requirement is 
found to be 



(<™> e 



H r Xf o T 



ITldm fcdm (Tfo ) ^dn 

7.6 lO^WV 1 - 



0.25 



'-H dm (Tfo) 20 



(57) 



This value is nearly a universal constant owing our good 
knowledge of the density parameter Sldm and the Hubble 
constant Hq. In particular it is almost independent of the 
Dark Matter particles' mass. Any Dark Matter particle 
that is able to annihilate must satisfy this condition in 
order not to overdose the Universe. This is true, in par- 
ticular, for any particle having a mass above the MeV 
range. 

The condition l|57(l however may be alleviated, to allow 
larger (but not smaller) values of (cv) . This is the case 
provided the particles and anti particles (with respective 
densities nd m + and nd m -) exhibits an asymmetry before 
the non-relativistic transition : 



ridr, 



ndn 



(58) 



Indeed, the density at the relativistic side then is nd m \r = 
ndm++ndm- and particle number conservation then shows 
the number density at the non-relativistic side to be 
n-dm\nr = idm+ — nd m -- To produce the observed relic 
number density, it may be shown (see Appendix IB. 3|) that 
a jump " dm l" r = a " r is required at the non-relativistic 

J i n dm \ r a cq{dm) 

transition. The relic density condition i|57|) then is replaced 
by the condition 



Q , eq(dm) 



(59) 



We now estimate the damping scale associated to Dark 
Matter candidates belonging to region II. At total decou- 
pling, Dark Matter particles have a velocity 



Vdm 



where 



^dec(dm) 



f 



Kdm(Tnr) 
Kdm (Tdec(dm) ) 

We then get 



l {II) /rd 

l sd I'drt 



vt 

7f \dec(dm) 



(60) 



(61) 



(62) 



The values we mention here slightly differ in case of co- 
annihilation. 



The ratio Pdm(Tdec{dm)) I ' Pr(Tdec(dm)) ~ a dec(dyn) I a eq IS 

still smaller than unity. But Td m -which depends on the ra- 
tio Pdm(Tdec{dm)) I ' $( T dec(dm)) of the species that are actu- 
ally coupled- may be unity if the Dark Matter is not cou- 
pled to radiation, as already noted for Region I. The self- 
damping and free-streaming damping lengths in Region II 
finally take the values given in Table 



Table 2. 



ii.ii) 



330 kpc rdm(T dec ( dm )) f g , 4 (Tdec(dm)) 



330 kpc f g', 3 {T dec(dm) ) 



1 dec(dm) 



The condition IjP < l s truct implies 



(T nr ) > 11 MeV/ 2 g'^(T dec{dm) ) 



dec(dm) 



{Is 



6 lCT^s- 1 \l00kpc) 



(63) 



4.2.3. Region III : a nr < a eq < a dec(dm) . 

As in the previous case, we have to consider the URFO 

and NRFO scenarios. To get the observed Dark Matter 

energy density one has to impose the same relations (|48|) . 

and (|57|l or ||5§J|, respectively, as for Region II. 

The velocity of the Dark Matter particles at decoupling 

has the same expression as for region II. A straightforward 

calculation (which takes into account that the decoupling 

occurs while the universe is matter-dominated) provides 

the damping scales. The interesting point here is that, 

i 

unlikely to region II, we expect the ratio Tdm ~ (Pdm/ f) 2 
to be of the order of unity (if not exactly equal to 1) 
since the density ^ refers only to the species coupled to 
the Dark Matter and the decoupling occurs in the matter 
dominated era. Therefore self-damping and free-streaming 
are comparable in this region: 



(in) 

fs 



— l sd I 1 dn 



~\dec(dm) i 



(64) 



which numerically yield the results given in Table [3J 
Despite / f s and l s d are given by the same expression than 
in Region II, their evaluation yields a different result since 
after equality, in Region III, the time of decoupling grows 
less rapidly with the interaction rate. As a result, the 
damping lengths turn out to be independent of this rate. 
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Table 3. 



Aim 



i(ni) 



435 



435 



kpc (- 
kpc (- 



l dm c{dm) 



)/l 



,(r nr ) ' 



19MeV 



0.3 



_1 // \- 2 

r2 ( L struct \ 

V lOOkpc J 



4.2.4. Region IV : a dec{dm) < a eq < a r , 



The Dark Matter particles are relativistic when they de- 
-i couple. The analytic expression of the damping scales 
is the same as in Region I. The free-streaming length 



(65) 



(IV) 



This provides a range for the Dark Matter mass : 

IH> draft dm (T nr ) > 



if g is still given by (|50|l . The numerical value of the 
collisional damping scale however is different. Since we 
have here pdm ~ the collisional damping length is still 
=smaller than the free-streaming length, but by a factor 

^\dec(dm) I ^\nr = 0>dec{dm) I a nr Only. 

Requiring the free-streaming damping scale to be 
smaller than I struct ~ lOOkpc still calls for Dark Matter 
particle masses in the keV range or above. This cannot be 
achieved in Region IV. Region IV is therefore excluded by 
the damping requirements. 



A specific case: coupling up to the onset of structure forma- 
tion. The previous expressions are only relevant provided 
decoupling occurs before structure formation, that is be- 
fore the non-linear collapse of the primordial structures, 
assumed to take place at a scale- factor a n i. We leave the 
latter arbitrary, but it is expected to be roughly of the 
order of ~ 1/10 for objects with a size < 1O 8 M . 

For a^ecUm) > a nh which may be rewritten as 
^dm\a=a n i > H„i (the interaction rate H n i is given in 
Table EH , there is no free-streaming at all before struc- 
ture formation. Only collisional damping is at work up to 
the non-linear collapse, with 



Isd/ I"dm 

that is 



(///') 



vr (H nl /\ 



435 kpc 



1 dm) — I 



(66) 



£2™ h 7 



0.3 

/ W < : - dm f ^ a r ) 



^dm{ r ^'dec(dm)) f 



lMeV 



- dm 



(67) 



It is useful to note that ij,^ 1 depends on the collision rate 
Tdm{T n i) taken at the time of the onset of the non- linear 
gravitational collapse and no longer at decoupling. Also, 

I 111 ' - 



lid 1 j ea.l|6T|l. is smaller than l 1 // , eq.Q 



although the corresponding Dark Matter interaction rate, 
taken at a given epoch, is larger. 

For the self-damping to be acceptable, we must have 

TUdmKdm 

(T nr ) > 19 MeV I m 70 ) r 2 d7n (T dec{drn ))f 2 



dm 



H 



nl 



( ^struct \ 

V lOOkpc J 



(68) 



Pending a more complete discussion of this issue in 
Paper II, we may recall that the present requirements are 
necessary conditions, as are all the bonds we establish in 
this paper, irrespectively of other astrophysical conditions 
which may be needed to achieve an acceptable scenario 
with the assumed Dark Matter parameters. 



4.2.5. Region V : a eq < a dec{dm) < a nr . 

The calculation is strictly identical to the one done in 
the previous section for Region IV since adec(dm) — a eq, 
does no longer correspond to a change in the expansion 
regime. The transition from a radiation to a matter dom- 
inated universe at Dark Matter decoupling occurs indeed 

for ddec(dm) — dnr- 

Region V is therefore also excluded by the damping 
requirements. 

4.2.6. Region VI : a eq < a nr < a dec{dm) 

Dark Matter decouples in the matter dominated era. The 
collisional, as well as the free-streaming scales are domi- 
nated by the contribution near the decoupling time. Since 
we have still p dm ~ their contribution is nearly the same 
as long as decoupling occurs before the epoch of non-linear 
gravitational collapse where the structures actually start 
to build up. The same calculation as in Region III is rele- 
vant in the present case. The conditions for the dm density 
to match the observed one still imposes (I48|l : anticipat- 
ing the astrophysical discussion of Paper II, we readily see 
that the URFO scenario is the only sensible case to con- 
sider. Indeed, in Region VI, the present-day number den- 
sity of these light Dark Matter particles is much higher 
than the photon density. Such an unusually large number 
density can only be reached under very special circum- 
stances. It does not make much sense to assume that it 
results from an even much larger one, which is followed by 
a period of annihilation. However, strictly speaking, this 
is not forbidden. 

The damping scales are given by their expressions in 
Region III, with comparable effects from collisions and 
free-streaming. Masses above the MeV scale are therefore 
required to avoid prohibitive damping. But such masses 
are not allowed in Region VI (fig-CJ. 

A specific case: coupling up to the onset of non-linear col- 
lapse. For r(T„;) > H n i, which we call Region VI', there 
is, on the other hand, never decoupling until the onset 
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of the non-linear gravitational collapse. Only collisional 
damping is at work so the free-streaming constraint drops 
out and expression (|67fl of the self-damping length, as well 
as the constraint (|68[) hold in this case. 

This leaves an allowed window in the parameter space 
of Region VP, for extremely large an interaction rate. 

5. Constraining the Dark Matter candidates from 
the neutrino induced-damping. 

The expression of the induced-damping scale shows 
that relativistic particles yield a large damping effect. 
They are expected to provide the most stringent con- 
straints on the Dark Matter interaction rates. Since the 
obvious components of the radiation are photons and neu- 
trinos, we specifically focus on these two species. These 
particles have the advantage that their dominant interac- 
tion (apart from the interaction with Dark Matter that we 
aim to discuss in the present paper) is with electrons. It 
is known and provides the scales of reference for the dis- 
cussion of the damping effect induced by these particles. 

In this section, we compute the collisional damping of 
Dark Matter fluctuations due to a possible coupling of 
Dark Matter with neutrinos and evaluate the correspond- 
ing constraints. The photon induced-damping of the Dark 
Matter fluctuations will be discussed in sectional 
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5.1.1. Dark Matter decouples from neutrinos while 
the latter are collisional. 

This is the case where 

^ dm — 1> \dec(dm—v) ■ 

The Dark Matter then experiences a collisional damping 
induced by its coupling to neutrinos. Since the neutrino 
interaction rate can be dominated by the collision rate 
with either the electrons or with Dark Matter, we are led 
to distinguish two sub-cases A and B. 

Region A. In this region, the neutrino interaction rate 
is dominated by interactions with electrons. This corre- 
sponds to : 

^u—dm — ^is — e\dec(dm — is) j that is ~ T^—g . (69) 

In this case, the neutrino collisional damping is due to 
their interaction with electrons. 

Region B. Here the neutrino interaction rate is dominated 
by interactions with Dark Matter. This corresponds to : 

dm ^ — e \dec(dm— u) t that is ~ ^v— dm • (^0) 

The collisional damping is then due to the neutrino inter- 
actions with Dark Matter. 



5.1. Dark Matter parameter space for neutrino 
induced damping. 

The neutrino induced-damping scale depends on the to- 
tal neutrino interaction rate with the medium rV = 
J2j IV— j = I^-e + T v -dm + In the standard scheme, 
this rate is dominated by the collisions with electrons, 

~ r„_ e . The only possibility for T v not to be dom- 
inated by T v - e is when T v -dm > r„_ e . 

Thanks to our approximations l|33[) . or (|41|l . the only 
other parameter which enters the calculation of the damp- 
ing scale is the upper limit of the integral l(39|l . namely 
the time td ec (dm-u) where Dark Matter decouples from 
neutrinos. This defines the reference time at which all our 
constraints on the interaction rate or on the cross-sections 
are obtained. 

For the actual applications, it turns out to be conve- 
nient to note that both the interaction rate T v ^dm and 
time tdec(dm-i/) niay be expressed as a function of the 
"reduced" (see Section [A. 2.6(1 interaction rate Tdm-v = 
Tdm-^a 3 , or equivalently as a function of {cfv) l/ _ dm — 
^v-dm/ridm- The other parameter entering the calculation 
is the Dark Matter particles' mass rridm- This defines a 
two-parameter space which allows one to classify all kinds 
of Dark Matter particles which interact with neutrinos, 
say [m dm ,Tdm-u] or [rridm, (ot)^^]. 

In the following, we shall distinguish "Regions" in this 
parameter-space which correspond to different analytical 
expressions of the neutrino induced-damping scale. All to- 
gether, they cover the whole parameter space. 



5.1.2. Dark Matter decouples from neutrinos while 
the latter are freely-propagating. 

In this situation, there is only one region, namely: 
Region C 

^dm—u ^ |c£ec(i/) ■ 

Neutrinos are already free-streaming when Dark Matter 
decouples from the latter. This is the mixed-damping 
regime. The corresponding damping length has been esti- 
mated ea. (|41|l . Since the neutrino collision rate does not 
enter the mixed-damping expression there is nq^point in 
the present case to separate the case where T v -dm is 
smaller or larger than r„_ e . 

5.1.3. Borderlines. 

The borderlines between regions A, B and C are shown 
in Figs. |2 E] El a nd depending on whether one con- 
siders the URFO or the NRFO scenario, and whether 
one displays the parameter-space in the [rridm, Tdm-v] ° r 
l m dm, (W)„_j m ] coordinates. These borderlines are calcu- 
lated explicitly in Appendix lC.2l They are important since 
they define the conditions of validity of our expressions of 
the damping length. 
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Fig. 2. Bounds in the [Dark Matter particles' mass, 
Dark Matter-neutrino interaction rate] parameter space 
obtained from neutrino induced-damping in the URFO 
scenario. The Regions A, B and C are distinguished by 
different colors. They correspond to different expressions, 
given in the text, of the neutrino contribution to the Dark 
Matter damping length. The dot-dashed lines separate the 
domains according to the ordering of the epoch of the Dark 
Matter - neutrino decoupling, the non-relativistic transi- 
tion and equality of the energy-densities. The hatches in- 
dicate in which part of the parameter space the neutrino 
induced-damping scale is greater than lOOfcpc (~ 1O 8 M0) 
scale. The indications in this figure are schematic. Some 
factors in general of order unity are omitted for simplicity. 
The reader interested by these constraints should use the 
expressions in the text, where all factors are given explic- 
itly. 
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Fig. 3. Bounds in the [Dark Matter particles' mass, Dark 
Matter-neutrino cross-section] parameter space obtained 
from neutrino induced-damping in the URFO scenario. 
The Regions A, B and C are distinguished by different col- 
orings. They correspond to different expressions, given in 
the text, of the neutrino contribution to the Dark Matter 
damping length. The dot-dashed lines separate the do- 
mains in the parameter space where the ordering of the 
Dark Matter - neutrino decoupling, the non-relativistic 
transition, the epoch of equality of the energy-densities 
changes. The hatches indicate the region in parameter 
space which is forbidden because the neutrino induced- 
damping yields damping above lOOkpc (~ 1O 8 M0) scale. 
The indications in this figure are schematic. Some factors 
in general of order unity are omitted for simplicity. The 
reader interested by these constraints should use the ex- 
pressions in the text, where all factors are given explicitly. 



5.2. Neutrino induced-damping scales and limits on 
the interaction rates. 

In the absence of significant Dark Matter interactions, as 
it is commonly assumed, neutrinos are expected to de- 
couple from the thermal bath in the radiation dominated 
era. The last interactions are with electrons. If we adopt 
the standard neutrino-electron interaction (see Section 
IA.3.1(I . the neutrinos decouple at a temperature slightly 
above lMeV, when their interaction rate is of the order 
°f r dec („_ e ) ~ 2.8 10 _30 s~ 1 . At neutrino decoupling, the 

damping scale is ldec(v-e) — Kct r r v de ^~ e > ~ 97 pc. 
The original calculation of this damping is due to Misner 
(1967). It corresponds to a mass scale somewhat below 
1 Mq. In the present case, the neutrino induced-damping 
scale takes this values in case the Dark Matter turns out 



to decouple from neutrinos just at the epoch the latter 
decouple from electrons. 

Let us now give the value of the damping lengths for 
all regions of the neutrino parameter-space. Analytical ex- 
pressions of these lengths may be found in Appendix ID. 21 
Pending the astrophysical discussion in Paper II, we may 
already notice that the case //> dm > $ makes little sense in 
case the relevant times to consider are close to the epoch 
of primordial nucleosynthesis, and is to be excluded. This 
is nevertheless possible somewhat before (provided annihi- 
lation occurs in time) or somewhat after (if Dark Matter 
particles are created late enough). 
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Fig. 4. Bounds in the [Dark Matter particles' mass, 
Dark Matter-neutrino interaction rate] parameter space 
obtained from neutrino induced-damping in the NRFO 
scenario. The Regions A, B and C are distinguished by dif- 
ferent colorings. They correspond to different expressions, 
given in the text, of the neutrino contribution to the Dark 
Matter damping length. The dot-dashed lines separate the 
domains in the parameter space where the ordering of the 
Dark Matter - neutrino decoupling, the non-relativistic 
transition, the epoch of equality of the energy-densities 
changes. The hatches indicate the region in parameter 
space which is forbidden because the neutrino induced- 
damping yields damping above the lOOfcpc (~ 10 8 M©) 
scale. The indications in this figure are schematic. Some 
factors in general of order unity are omitted for simplicity. 
The reader interested by these constraints should use the 
expressions in the text, where all factors are given explic- 
itly. 
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Fig. 5. Bounds in the [Dark Matter particles' mass, Dark 
Matter-neutrino cross-section] parameter space obtained 
from neutrino induced-damping in the NRFO scenario. 
The Regions A, B and C are distinguished by different col- 
orings. They correspond to different expressions, given in 
the text, of the neutrino contribution to the Dark Matter 
damping length. The dot-dashed lines separate the do- 
mains in the parameter space where the ordering of the 
Dark Matter - neutrino decoupling, the non-relativistic 
transition, the epoch of equality of the energy-densities 
changes. The hatches indicate the region in parameter 
space which is forbidden because the neutrino induced- 
damping yields damping above the lOOfcpc (~ 1O 8 M0) 
scale. The indications in this figure are schematic. Some 
factors in general of order unity are omitted for simplicity. 
The reader interested by these constraints should use the 
expressions in the text, where all factors are given explic- 
itly. 



5.2.1. Expressions of the neutrino-induced damping 
scales in Region A and B. 



Region A : 

t / dm— v \ ct . 

V r„_ e ) a 

or more explicitly 



lud = 97 pc r v ng 



2.8 10" 30 s 



30 o -l 



(71) 



(72) 



The largest value that this damping length can take is 
when td ec {dm-u) = tdec(v-e)- This case only provides rel- 
evant limits if we require primordial fluctuations to exist 



down to very small scales of less than 100 pc (~ 0.1 M©). 
It can be obtained from ea. (|72j) by writing l v d < l s truct 
but will not be evaluated explicitly here as not of direct 
cosmological interest. 



Region B 



lvd = 



dva—u 
u—dm 



ct 

^ \dec(dm — u) 

a 



(73) 



Only the case of a radiation-dominated universe is rele- 
vant here since the converse leads to prohibitive damp- 
ing. Hence, as is readily seen in figs. to (J5J, we need 
to consider only the case where Dark Matter is rela- 
tivistic {adec(dm-v) < o-nr)- With the symmetry relations 
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of the interaction rates (Sections 12.41 and lrj.2f> . we get 
the explicit relations displayed in Table 0] and [SJ Note 
that we do not consider explicitly the NRFO scenario for 
a-eq(dm-v) < 0-nr, a case which is possible but quite un- 
likely for the reasons discussed above, see Section 14.2.61 

Also, for the neutrino interaction rates relevant to our 
calculation, decoupling of the Dark Matter with neutri- 
nos is to occur well before the standard epoch of matter- 
radiation equality and the fate of the neutrinos afterwards 
is irrelevant to our purpose: with the present limits on neu- 
trino masses we remain well within the epoch where the 
neutrinos are fully relativistic. 



again be transformed into limits for {o'v) 1J _ dm by means 
of the symmetry relations between the interaction rates 
(Sections E3 and BO)) . 

This calculation does not require any knowledge of the 
number of electrons which are present, and holds whether 
the latter have annihilated or not. Also, all times of con- 
cern are well before the standard epoch of mater-radiation 
equality. It is worth to remember at this stage that our 
calculation are valid for massive neutrinos, provided their 
mass is within the known physical and astrophysical lim- 
its. 



Table 4. Region B (URFO scenario) 



Table 6. Region C (URFO scenario) 



'7» <v(-^) V. ; (%^)" 

«dm \ 2 T dm-v 

eV J 2.8 10-™ s- 1 
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( "dmfeTO 2 A 2 ( m dm K, dm \- 
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Table 5. Region B (NRFO scenario) 



97 pc r 



2.8 10~ i0 s- 



Table 7. Region C (NRFO scenario) 



-1 -1 ( 9, dm \ 2 I l.tmt 
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Ivd 
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lOOfcpc 



< 9.0 10" 28 cm 3 s _1 



r -2 I 4g' ,„ \ „l\ m, dm n dm I l struat V 
v \ 3 J y * lMeV V lOOfcpc J 



5.2.2. Expressions of the neutrino-induced damping 
scales in Region C. 

The mixed-damping regime is at work in this region. We 
have (EqEJ 



ct 



Ivd ~ TtTj/ \dec(dm— v) 



(74) 



Limits on Tdm-u are given in table and respectively, 
for the URFO and the NRFO scenario. These limits can 



^dec(dm-i') ^ a nr 

does not exist in this case 
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Fig. 6. Bounds in the [Dark Matter particles' mass, Dark 
Matter -photon interaction rate] parameter space obtained 
from photon induced-damping in the URFO scenario. The 
Regions A and B are distinguished by different colorings. 
They correspond to different expressions, given in the text, 
of the photon contribution to the Dark Matter damp- 
ing length. The dot-dashed lines separate the domains 
in the parameter space where the ordering of the Dark 
Matter - photon decoupling, the non-relativistic transi- 
tion, the epoch of equality of the energy-densities changes. 
The hatches indicate the region in parameter space which 
is forbidden because the photon induced-damping yields 
damping above the lOOfcpc (~ 1O 8 M0) scale. The indica- 
tions in this figure are schematic. Some factors in general 
of order unity are omitted for simplicity. The reader in- 
terested by these constraints should use the expressions in 
the text, where all factors are given explicitly. 
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Fig. 7. Bounds in the [Dark particles' Matter mass, Dark 
Matter -photon cross-section] parameter space obtained 
from photon induced-damping in the URFO scenario. The 
Regions A and B are distinguished by different colorings. 
They correspond to different expressions, given in the text, 
of the photon contribution to the Dark Matter damp- 
ing length. The dot-dashed lines separate the domains 
in the parameter space where the ordering of the Dark 
Matter - photon decoupling, the non-relativistic transi- 
tion, the epoch of equality of the energy-densities changes. 
The hatches indicate the region in parameter space which 
is forbidden because the photon induced-damping yields 
damping above the lOOfcpc (~ 10 8 M©) scale. The indica- 
tions in this figure are schematic. Some factors in general 
of order unity are omitted for simplicity. The reader in- 
terested by these constraints should use the expressions in 
the text, where all factors are given explicitly. 



6. Constraining Dark Matter properties from the 
photon induced-damping scale. 

As already mentioned, the possible coupling between Dark 
Matter and photons also induces a source of collisional 
damping. The associated damping scale can be in- 
ferred from the formula (|33|l . The calculation, however, 
is somewhat different from the one for l v d because the 
epoch at which the Universe becomes matter dominated 
and the epoch of the recombination now get into play. 
This introduces two new scales H eq ^ m ^^ and H rec which 
are the "reduced" (H = Ha 3 ) Hubble rates at some 
specific time related to the epoch of equality and re- 
combination, respectively (precise definition are given in 
Sect. IA.3.3(l . Similarly to the neutrino case, from T 7 = 
d m + .... we see that the relevant interaction 



rates are T 7 , T 7 _ e and r 7 _d m (see Section lA.2.6|) . So, the 
parameter space for the photon induced-damping may be 
taken as [m dm , Fi m _ 7 ] or [mdm, fiv) 1 - dm ]- 

6.1. Dark Matter parameter space 

"Regions" are expected to appear, as in the neutrino case, 
depending on whether photons are collisional or collision- 
less. The new discussion here is about the role of recom- 
bination and the relevance of the epoch of matter dom- 
inance. The borderlines of these regions are given ex- 
plicitely in Appendix IC. 31 
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Fig. 8. Bounds in the [Dark Matter paticle's mass, Dark 
Matter -photon interaction rate] parameter space obtained 
from photon induced-damping in the NRFO scenario. The 
Regions A and B are distinguished by different colorings. 
They correspond to different expressions, given in the text, 
of the photon contribution to the Dark Matter damp- 
ing length. The dot-dashed lines separate the domains 
in the parameter space where the ordering of the Dark 
Matter - photon decoupling, the non-relativistic transi- 
tion, the epoch of equality of the energy-densities changes. 
The hatches indicate the region in parameter space which 
is forbidden because the photon induced-damping yields 
damping above the lOOfcpc (~ 1O 8 M0) scale. The indica- 
tions in this figure are schematic. Some factors in general 
of order unity are omitted for simplicity. The reader in- 
terested by these constraints should use the expressions in 
the text, where all factors are given explicitly. 
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Fig. 9. Bounds in the [Dark Matter particles' mass, Dark 
Matter -photon cross-section] parameter space obtained 
from photon induced-damping in the NRFO scenario. The 
Regions A and B are distinguished by different colorings. 
They correspond to different expressions, given in the text, 
of the photon contribution to the Dark Matter damp- 
ing length. The dot-dashed lines separate the domains 
in the parameter space where the ordering of the Dark 
Matter - photon decoupling, the non-relativistic transi- 
tion, the epoch of equality of the energy-densities changes. 
The hatches indicate the region in parameter space which 
is forbidden because the photon induced-damping yields 
damping above the lOOfcpc (~ 10 8 M©) scale. The indica- 
tions in this figure are schematic. Some factors in general 
of order unity are omitted for simplicity. The reader in- 
terested by these constraints should use the expressions in 
the text, where all factors are given explicitly. 



6.1.1. Dark Matter decouples from photons while the 
latter are collisional. 

This implies 

Fdm — 7 5; 1^7 \dec(dm— 7) • 

Region A. Photons collisions dominated by photon- 
electron scattering. 

1^7— dm — ^7— e \dec(dm— 7) 

Region B. Photons collisions dominated by photon- Dark 
Matter scattering. 

^7— dm !\ — e \dec(dm— 7) 

The separation between the two regions A and B may, 
among others, occur for dm — 7 decoupling before recom- 



bination, that is for a sizeable T 7 _ e rate. This requires 
simultaneously sufficiently small Tdm-'y rate for the de- 
coupling be early enough, and a sufficiently large r 7 _d m 
rate to be larger than the 7 — e interaction. These two, 
somewhat antagonistic, requirements can actually be sat- 
isfied (see the evaluation in Appendix IC.3(I . but for unre- 
alistic values of the Dark Matter particles' mass mdm (of 
the order of 10~ 2 eV). This separation is barely visible in 
figs. ® and (0. 

More realistically, the separation between the regions 
A and B simply occurs at recombination: before the latter 
the 7 — e interactions dominate, and after, only 7 — dm 
interactions are left. So, Regions A and B in practical 
cases correspond just to the separation between dm — 7 
decoupling before and after recombination (Table 01. 
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6.1.2. Dark Matter decouples from photons while the 
latter are freely-propagating. 

This is the analogue of region C in the neutrino case. It 
requires r 7 < H at the Dark Matter -photon decoupling 
(Tdm--y = H). The earliest this may happen, whatever 
the strength of the Dark Matter interactions, is at re- 
combination: so, Dark Matter -photon decoupling neces- 
sarily occurs at an epoch where the universe is matter- 
dominated. One then has necessarily Td m --y < I^-dro ac- 
cording to the relation l|21|) . But the first relation implies 
I^y-dm < H and in turn the second Tdm-^ < H. Therefore 
Dark Matter cannot remain coupled to photons while the 
latter are freely propagating. So there is no Region C, i.e. 
no mixed-damping regime. 



We get in this case 

i 7 d = 7Tr 7 ( ~" 7 ) — \dec(dm-y) ■ (75) 

V r Th J a 

We are led to consider the two following cases : 

1. Dark Matter decouples from photons in a radiation 
dominated Universe, that is r<j m _ 7 < H eq (d m -~f)- In 
this case the shear viscosity dominates. This yields 
the damping length as well as the constraints given 
in Table El 

Table 9. Region A. Decoupling in a radiation dominated 

Universe (a dec(dm _ 7 ) < a eq ) 



Table 8. Definition of the two Dark Matter Regions useful 
for the photon induced-damping calculation 













Region A 


^ dm — 7 


^ T-y — e \dec{dm — 7) 


^dec(dm- 


-7) — &rec 


Region B 




~22i dm \dec{dm — 7) 


^dec(dm- 


-7) ^ Oir&c 



6.2. Photon induced-damping scales and limits on the 
interaction rates. 

We now compute the damping scales associated with the 
regions A and B for the two scenarios NRFO and URFO. 
Region A is slightly more complicated than region B be- 
cause it can potentially exist in both a radiation and mat- 
ter dominated Universe while region B necessarily takes 
place in a matter dominated Universe. We give here the 
most relevant numerical results. A more systematic pre- 
sentation of our analytical results is given in Appendix 

EH 

6.2.1. Expression of the photon-induced damping scale 
for Region A. 

The condition T 7 _ e > r 7 _<j m is in some sense the stan- 
dard situation since, in this case, the photon decoupling is 
set by the recombination epoch. The damping of the Dark 
Matter fluctuations is then induced by the interaction of 
photons with the electrons. Below a temperature around 

100 keV (a case which is relevant for Td m f > 10~ 28 s _1 

and is the only regime of practical interest), the photon- 
electron elastic scattering cross-section (i.e. the Thomson 
cross-section) yields the reduced interaction rate T Th = 
er Th c n e ^5 10~ 21 s _1 , while the recombination cor- 
responds to a (reduced) Hubble rate H rec ~ 3 10 _23 s _1 



1 / \ 

1,4 = 8.2 MpcrW? ( 61 r o d -"«7-i j 
f dm _ 7 < 3.2 10- 25 s" 1 r~i g'l (^)'(^) 1 

&dec(dm — 7) ^* Q nr 

«dm<^J) 7 _ dm < 1.1 10" 23 cm 3 8- 1 r~tg'${T) 

( n b h 70 2 \3 t m dm K dm \ I Struct \ 3 
I 0.05 j \ lMeV J \100kpcJ 

a dec(dm — ~() < a nr 

Kdm(av)^ dm < 5.1 lO" 28 em's- 1 r^g'l (^f 2 -) 

( (l b h 7a 2 \ 3 ( l st ruct \l 
I 0.05 J V lOOfcpc ) 



2. Dark Matter decouples from photons^ in a matter dom- 
inated Universe, that is Tdm-y > H eq (d m -i)- This is 
the case where the heat conduction dominates, so that 

V = 8.2M P cr 7 (^) I (W)" 1 5 (^)" 1 




This scale is typically above IMpc. So, when the 
Dark Matter -photon decoupling is in the matter 
dominated era, the photon induced-damping is always 
prohibitive according to our criteria. 

It may be of interest to note that for the largest 
interaction rate allowed in Region A, that is for 
r<im-7 = H rec , i.e. the Dark Matter decouples from 
the photons at the recombination epoch. The scale l^d 

thus reaches Is = r Kct m r 1 (^H rec a r ec/^Thj , the Silk 
damping length which may be estimated as 

l s = 4.9 Mpc (Q m h 70 2 )-i (n b h 70 2 )'i . (77) 
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This length corresponds to a total mass (including the 
Dark Matter ) 

M s = 6.6 10 13 Mq (n m h m 2 )i{n b h 70 2 y? . (78) 

It may be noted that, with the now well established 
values of fi m , Ho and especially Qb, these are very large 
scales, of the order of Is — 30 Mpc and Ms = 10 15 M©. 



6.2.2. Expression of the photon-induced damping scale 
for Region B. 

The case where Dark Matter decouples from the pho- 
tons after the recombination obviously leads to a damping 
which is larger than the Silk damping. We shortly treat 
this case for completeness. For cross-sections just above 
this threshold, that is for Tdm-y ~ H rec , the reduced 
photon- Dark Matter interaction rate is required to be 
r 7 -dm = -ff e9 (dm-7)(-ffrec/-ff eg (dm- 7 )) 5/3 and is smaller 
than T Th . Since, after recombination, the photon-electron 
interactions become negligible, only T^-dm contributes to 



T 7 . Such a photon- Dark Matter interaction therefore in- 
duces a sudden increase of the transparency (and of the 
viscosity coefficients), inducing in turn an increase of the 
damping length : 



l^d 



'v V 2 * 

-L rim — 7 \ . 

^ ■y — dm f 



(79) 



More specifically, we obtain (for non-relativistic Dark 
Matter, as is the case in nearly all of Region B): 



jd 



350 Mpc r 7 



0.3 



0.25 



(80) 



The damping length stays constant (and equal to the 
above value) whatever the interaction rate. For larger 
dm — 7 interaction rates, the viscosity coefficient gets 
smaller but the interaction occurs a longer time. The two 
effects happen to exactly compensate. 

The damping in this case, obviously, is prohibitively 
large. This is also expected to be the case for the vanish- 
ingly small region in the parameter-space where the Dark 
Matter is relativistic. 

6.3. Limits on indirect Dark Matter -photon 
interactions. 

Even neutral Dark Matter may experience photon 
induced-damping. Indeed, if the Dark Matter happens 
to be directly coupled to any of the species j which are 
themselves coupled to the photon fluid (baryons, electrons 
etc.), the above damping is relevant. Photon induced- 
damping then is at work as long as the Dark Matter is 
coupled to the fluid j, and the previous estimates are rel- 
evant, with the replacement 



This can be used to get constraints on the relevant dm — j 
cross-sections. We for instance get 



(™)dm-b < 1-6 10~ 18 cmV 1 r T fy£(T) 



flhh 



bn-70 



2\ S 



(Istruct \ 3 
lOOfcpc J 



(82) 



dm— 7 



dm—j 



(<rv) 



dm-j n j 



(81) 



Exactly the same constraint holds for the dm — e cross- 
section. 

7. Conclusions. 

In this general evaluation of damping properties, we have 
developed an approach which allows to discuss systemat- 
ically all possible kinds of Dark Matter. 

We have written the transport coefficients for a com- 
posite fluid of species in terms of the interaction rate of 
each of these species with the medium. This is rather 
straightforward for the shear viscosity and can be read- 
ily extended to heat conduction. We have taken some care 
to discuss bulk viscosity along the same lines, although 
the latter turns out for the sake of the present paper not 
to generate stronger constraints on the collision rates and 
cross-sections than shear viscosity. This provides a new 
insight into the physical processes at work, and has been 
used as an opportunity to correct some inaccuracies and 
wrong statements in earlier papers. We finally wrote all 
the transport coefficients of a composite fluid as a sum 
of contributions from the individual components, taking 
the coupling to the whole medium into account. This is 
the key step which allows us to discuss the damping of all 
possible kinds of Dark Matter. 

We are led to separate, following iBoehm et all (20011. 
the contribution due to self-damping, that is the transport 
by the Dark Matter itself, and the one due to the induced- 
damping which is the transport by other species effective 
for damping the Dark Matter fluctuations. The classical 
examples of which being neutrinos and photons. This work 
aims to put our previous results on a general and firm 
basis. 

Writing that each of these contributions alone must 
not produce prohibitive damping, we were able to get 
necessary conditions separately for the Dark Matter in- 
teraction rate with all coupled species, as well as for its 
specific interaction with each of the species it is coupled 
to. More generally, we identify the parameter space of rel- 
evant, namely interaction rate and Dark Matter particles' 
mass. Another parameter, which is seen to appear at the 
proper places, comes into play. The latter (K c i m in the 
text) measures the possible offset of the Dark Matter and 
photon temperatures. It is in practice of order unity or 
somewhat larger, but it must be stressed that in special 
circumstances (e.g. axion-like Dark Matter) this parame- 
ter may get overwhelmingly large and provide the domi- 
nant physics. 

As expected, the bounds due to self-damping are 
closely related to the damping due to Dark Matter free- 
streaming. Both bear similar properties: the velocity of the 
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Dark Matter particles, that is their mass, is the most im- 
portant parameter. However, there are regions of parame- 
ter space where the Dark Matter interaction plays a role. 
As a result, in these regions, the expression of the free- 
streaming length (expressed in term of the Dark Matter 
mass), which is generally used for WDM candidates, does 
not hold anymore. Instead the free-streaming length de- 
pends on both the mass and interaction rate of the Dark 
Matter candidate. Said differently, the well-known require- 
ment of having Dark Matter particles with a mass above 
~ IkeV is valid only in a restricted region of parame- 
ter space, which covers only a fraction of the Dark Matter 
models, where the Dark Matter interactions are really neg- 
ligible. In regions where the Dark Matter interaction rate 
is not too large, the smallest acceptable value of the Dark 
Matter mass is above lMeV. In other parts, where the in- 
teraction rate turns out to be very large, extremely small 
masses (notably in the eV range) are allowed. 

The bound due to the damping induced by a given 
species merely depends on the interaction rate of the Dark 
Matter with this species, as well as on the species' particle 
velocity. Relativistic species are the more efficient, damp- 
ing the Dark Matter fluctuations nearly over the horizon 
size. 

The neutrino induced-damping yields quite a rich 
parameter-space diagram (Dark-Matter-neutrino interac- 
tion rate/ Dark Matter particles' mass). The bounds due 
to the possible coupling of the Dark Matter to the neu- 
trinos turn out to be quite novel. Indeed, one in partic- 
ular refers to Dark Matter that is coupled to neutrinos 
which are already free-streaming! We have evaluated an- 
alytically this effect in what we have called the mixed- 
damping. The difference with the collisional bound, which 
is only a lower bound to the damping length but takes only 
into account the effect of standard collisions, is large. Our 
mixed-damping estimate is still rather rough, although 
quite natural a result. It definitely would deserve a check 
by numerical simulations. The difference with the colli- 
sional bound is so large that there should be a substan- 
tial effect even in a realistic calculation, and the (lamping 
of a collisional fluid coupled to the free-streaming neu- 
trinos should be the dominant contribution as we have 
stressed here. This is important, especially in the context 
where Warm Dark Matter scenarios. Indeed, we saw that, 
when our constraints are marginally satisfied, the damp- 
ing induced by the neutrinos could yield a collisional Dark 
Matter candidate. 

The bound due to a possible coupling of Dark Matter 
with photons (this may be the case through higher or- 
der reactions involving charged particles even if the Dark 
Matter is neutral), in the (Dark-Matter-photon interac- 
tion rate/ Dark Matter particles' mass) parameter space, 
are somewhat less constraining. They are nevertheless rel- 
evant for a large fraction of the parameter space. Again, 
the possibility for Warm Dark Matter scenarios, where the 
damping is not due to free-streaming, but to interactions, 
exists. 



The major (and nearly the only) hypothesis is that 
the Dark Matter went through a phase of statistical equi- 
librium, at some stage of its evolution. There is also the 
(milder) assumption that a sudden decrease by many or- 
ders of magnitude in the Dark Matter particle density is 
not possible elsewhere than at its non-relativistic transi- 
tion (for instance a phenomenon like the recombination 
for the electrons -where the latter are substracted from 
the medium at a temperature bearing no relation with the 
electron mass- is assumed not to occur for Dark Matter ). 

Our bounds on the Dark Matter interactions with the 
most relevant species have been calculated without any 
assumption on the form of the primordial spectrum. The 
calculations are valid for any Dark Matter interaction rate 
and any Dark Matter particle mass, the values of practical 
interest covering 20 orders of magnitude for the interaction 
rates or cross-sections and much over 13 orders of magni- 
tude for the Dark Matter particles' mass. Our estimates 
are valid within factors of order unity, a procedure which 
provides important simplifications in the classification. A 
thorough discussion of all possibilities allows to place any 
scenario within the relevant parameter space, with a pre- 
cise statement on the important parameters and the ac- 
tual physical process at work. The reader interested in a 
specific scenario can easily work out the relevant bound 
with more accuracy. He/she will find here a useful esti- 
mate of the result and a guide on how to proceed. The 
present results provide, we believe, a new starting basis 
for a systematic discussion on the possible nature of Dark 
Matter. The astrophysical relevance of these findings, in- 
cluding other constraints than just the damping condition, 
or more restrictive as sumptions on the in itial conditions, 
in particular those of l|Boehm et alJ l)2002lU which specifi- 
cally restrict themselves to CDM initial conditions in order 
to yield more stringent limits, will be thoroughly discussed 
in Paper II. 
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Appendix A: Definitions and notations. 

We give here the expressions actually used in the text. 
They may be obtaine d from elementary textbooks (e.g. 
iKolb fc Tiiriierl \l<W(t) ). 



A.l. Present-day Universe. 

We write the present-day Hubble constant as 
Hq = 70 h 70 km s~ l Mpc^ 1 . 

A. 1.1. Temperature. 

We use for the present CMB temperature the value 
To = 2.73 K = 2.35 10~ 10 MeV , 
that is 

he 9 

— = 8.39 10~ 2 cm . 

Jo 

A. 1.2. Energy densities. 

The present critical, photon and neutrino (the latter for 
our "reference" cosmology) energy densities are 

p 0c = 9.21 10~ 30 h 70 2 g cm" 3 

= 5.16 10~ 3 h 70 2 MeV/c 2 cm" 3 
= 1.36 10 11 h 70 2 M G Mpc- 3 

= 4.67 10~ 34 g cm~ 3 
p 0v re f = 3.19 10" 34 g cm~ 3 
Po( 7 +u rof) = 7.86 10~ 34 g cm' 3 . 

The ratio of the total matter, Dark Matter and baryon 
energy densities to the critical density are denoted respec- 
tively Q m , fl dm and Q b : 

POm = ^mPOc 
POdm = QdmPOc 
POb — QbPOc ■ 

A. 1.3. Number densities. 

The Dark Matter , photon, neutrino (the latter in our 
"reference" cosmology), and baryon number densities are 
respectively equal to : 

no7 = 413 cm~ 3 
nou rcf = 338 cm -3 

i on in-3 -3^dmh 7<) ( md 

nodm — 1-29 10 cm ' 



( m dm \- 1 
\\Mev) 



nob = 2.75 10 7 cm 



-3 ^bh 7a 



0.25 \lMeV< 

2 



0.05 



The latter value has been obtained assuming an average 
mass per baryon of 938Af eV, that is the proton mass. It 
will also be used as the comoving number of baryons since 
we neglect in this paper the slight average mass offset, 
upwards before nucleosynthesis, and downwards after. 
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A. 2. Evolution. 
A. 2.1. Times. 

We write the scale factor as 

' t 

This expression holds even if is not unity, and/or A non- 
zero provided the scale factor a is not too close to unity, 
say for a < |, if fl is not much smaller than 0.2, and A not 
much larger than 1 . The parameter a is equal to 1 /2 in the 
radiation dominated era and equal to 2/3 in the matter 
dominated era. The normalization t a may be obtained by 
matching the above expression to the Friedman equations. 
In the radiation dominated era, t a is equal to 



t r — El/2 — 



16ttG ff »(T) 



P~f(T ) 



3 k 4 (T) 
2.32 10 5 Mpc/c g'~*{T) 



= 2.3910 19 s g\ 



where (see Sect lA.2l4l for more detail) g' t 
the matter dominated era, t a is equal to 



<?*/3.36. In 



*m=*2/3 = [67rG p m (T ) 



5.3710 17 a 



0.3 



5.21 10 3 Mpc/c 



0.3 



At the epoch of equality of matter and radiation, the 
expansion factor is not a power-law and the above approx- 
imations break down. They have been established so as to 
be exact far from this point. For this reason, it is slightly 
inaccurate to define the time of equality by equating the 
two expressions of the expansion parameter one obtains 
for a = 1/2 and a = 2/3. A better way of defining the 
epoch of equality in terms of the expansion factor a eq will 
be given in the following section, equating the associated 
Hubble rates which are more directly connected with the 
energy densities of the various species. With the definition 
of a eq given in Section lA.2.21 we have 
i 

A. 2. 2. Hubble rates. 

We introduce the constant 
a 



H, 



The latter is in the radiation dominated era equal to 

H r = Hx/2 = 2.10 10- 2( Vyj(T) 

and in the matter dominated era equal to 



H m = H 2/3 = 1.24 10- 18 s 



0.3 



This allows us to write the relation between the expansion 
factor and the Hubble parameter at a given epoch as 

H 

H a , 

It will also be convenient, for further use, to define a "re- 
duced" Hubble rate 

H = Ha 3 . 

Hence we have also 




The epoch of matter-radiation equality will be defined by 
equating the two expressions of a obtained for a = 1/2 
and a = 2/3, so that the Hubble rate, which is propor- 
tional via the Friedman equations to the energy-densities, 
is continuous at the transition, reflecting the continuity of 
the energy densities. We then get at the epoch of equality 

2v 2 



n n,,, = = 2.59io- 13 s-V; s (r eg ) 



and 



H 3 
eq ~ H 2 

rn 

This implies 
H, ^ 

or, equivalently 



H = H, 



§,miQ- M s- x g'l{T eq ) 



0.3 



0.3 



A. 2. 3. Temperatures. 

The temperature of a given species i with a thermal dis- 
tribution 9 will be denoted Tj. We omit the index i when 
dealing with photons. For each species i, we may define a 
parameter Ki (T) which relates the temperature Tj of the 
species i to the scale-factor a : 

T- = — 
Ki(T) a 

This parameter Ki is constant in most cases. Its value at a 
given time depends on the history of the system. It can be 
determined by using the comoving entropy conservation 
when the latter is indeed conserved. As is well-known, the 
term Ki changes if, for instance, the fluid to which i be- 
longs contains another species which is under way of an- 
nihilating. We use k 7 (T) = k(T), The latter is equal to 
unity nowadays : k(Tq) = 1, while k v (Tq), which is asso- 
ciated with the present neutrino temperature, is, in case 
our "reference" cosmology is adopted, equal to (11/4)3. 

9 Even if the distribution is non-thermal, one can define by 
averaging over the actual distribution an effective temperature 
Ti —< |j >, and accordingly a parameter «i(T), which for 
most if not all applications in the present paper, play the same 
role. 
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A. 2. 4. Energy densities. 

The energy density of a relativistic species is given by : 

9i rri 

Pi = £^a SB l l , 

with a SB is the Stephan-Boltzmann constant, gi the spin- 
degeneracy factor and e a factor equal to 7/8 or 1 for, 
respectively, Fermi or Bose-Einstein statistics of relativis- 
tic particles with negligible chemical potential. 10 

It turns out to be convenient to write the previous 
relation as 



3*1 rri4 

Pi = -^-as B J , 
where 



9*i = £9i — 



We also define the quantity 

g*i(T) 



3.36 k 4 (T) ' 



which is more useful than g„ for the calculations of the 
damping scales. This yields 



/ /jt>\ P0(j+u ref) 

Pi =9 *i( T )— -hi ■ 



The parameter g'^ is, under normal circumstances, of the 
order of unity. We have, at the present epoch for our "ref- 
erence" cosmology, g\ v rcf = 0.41 and g' „ 7 ref = 0.59. 

We also have, summing over all relativistic species in 
the universe 

5 '*cn = v 5 UT)= •"' ,Ji 



3.36 k 4 (T) 



so as to get 

I (rps P<d(l+v ref) 
Pr =9 A T )— ^-4 L ■ 

In our " reference" cosmology, this factor is unity after the 
electron annihilation. At early times, g'*(T) behaves as 
1.9 g: 1/3 (T). 

Finally, we define 

h = Pi+Pi 

where pi is the pressure of species i. For relativistic matter, 
p'i^'SPi while for non-relativistic matter p ! , i ~ pi. 



10 The factor e may differ from 7/8 or unity if we consider, for 
instance, relativistic particles with non-zero chemical potential 
(as may be the case if they are chemically decoupled but still 
in thermal equilibrium). 



A. 2. 5. Number density 

We use throughout this paper comoving number densities, 
namely : 

n = n a 3 . 

The latter appear more convenient because they are con- 
stant as a function of time when the particle number is 
conserved, so we keep only the dependence due to the mi- 
crophysics, the cosmological dependence being removed. 
In particular, the comoving number density of protons or 
electrons before or after their non relativistic transition, 
are respectively 

n p (T b > m.p/3) = n e (T e > m e /3) , 

n b (T b < m. p /3) — n p (T b < m p /2>) = n e (T e < m e /3) . 

For the electron density, we neglect the slight offset due to 
the presence of Helium (that is neutrons) after the epoch 
of nucleosynthesis. Clearly also, the electron density is rel- 
ative to free electrons only before recombination. 

For a relativistic species i, including possibly the Dark 
Matter itself, it will turn out to be more convenient to use 
(e.g. in Appendix IB. 21 to relate various interaction rates) 
the quantity 



3Tn 



625 cm 



-3 V„(T) 



rather than the comoving number density n\ of species i, 
to which however rf, i is close. 

A. 2. 6. Reduced interaction rates. 

As for the comoving number density, to remove the most 
obvious cosmological dependence of the interaction rates T 
and leave only the temperature dependence related to the 
microphysics, we introduce "reduced" interaction rates : 

f = T a 3 . 

We define the time t dec ( dm _ i - ) and the corresponding 
scale-factor a ( j ec ( ( j m _ i ), related to the epoch where Dark 
Matter particles decouple from the species i. It is, by con- 
vention, taken to be the epoch where the interaction rate 
of the Dark Matter with the species i, namely T^m-i is 
equal to the Hubble rate : 



dra—i 



= H 



The value of the reduced interaction rate T^m-i at the 
Dark Matter decoupling from the species i will be denoted 



it r 



dec(dm-i) 



It corresponds to 



ec(dm- i } 



— V 

— J- dm-i M \a=a dec{dm _ i) 



We also need the total interaction rate of a species i, 
which includes the interactions with all species j (includ- 
ing i) 
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The quantity T dec ^ will then denote the value of 1^ at the 

where species i totally 



time tdecU) and scale-factor a de 
decouples (even with itself) . 



For generality, we denote by T 



dec(x) 



and T 



dec(x) 



the 



reduced and normal interaction rate associated with x, 
where x stands either for a particular coupling i — j be- 
tween the species i and j or simply for the full coupling 
of a species i to the medium. The relation between T dec ^ 



and T 



dec(x) 



may be written 



^dec(x) H a 



- dec(x) 



This relation takes two different forms 
Radiation dominated era : 



r ' dec(x) — H e q 



H, 



Matter dominated era 



r ' dec(x) — H e q 



He, 



A. 3. Typical scales. 

A. 3.1. Interaction rates. 

The neutrino-electron and the photon-electron interaction 
rate play a central role in the determination of the damp- 
ing induced by respectively the neutrinos and the photons. 

Taking for the neutrino-electron cross-section (aver- 
aged over particles and antiparticles as well as over the 
three flavors) 



cr„_ e c = 4.02 KT^cmV 1 



EuE„ 



(IMeV) 2 
we get, after statistical average 



r — 



: v) n e 



where n e is the total (relativistic) electron density includ- 
ing electrons and antielectrons. This rate varies as oc T 2 
because of the energy-dependence of the cross-section. At 
v — e decoupling (given by r„_ e = H), this rate takes the 
value T dec{v _ e) . 

For non-relativistic electrons, up to recombination, the 
7 — e interaction rate is a constant as a function of time : 



r 

-L -y — e 



r. 



cr Th c n e 



where 



3m; 



is the Thomson cross-section. 



Numerical values of these rates are given in table IA.21 
A. 3. 2. Scale-factors fTablelA~Hl. 

Several typical scale-factors may be defined. They are as- 
sociated with the time-scales used for the calculation of 
the damping effects. 



There first is a nr , the epoch at which the Dark Matter 
becomes non-relativistic, taken by convention to corre- 
spond to T dm = m dm /3. 



We are led to define the scale-factor a, 



eq{^-\-v ref) 



which 



characterizes the epoch of equality between matter and 
photon+neutrino energy-densities in our "reference" cos- 
mology and a eq the true epoch of matter-radiation equal- 
ity in the general case, related to H r and H m already 
defined in Section lA. 2. 21 

In order to get simpler expressions for the borderlines 
in the induced-damping case, we define also 



eq(dn 



i) 



the ratio of the energy-density of a relativistic species i 
extrapolated to the present-day to the present-day energy- 
density of the Dark Matter . If the Dark Matter is non- 



relativistic at a, 



q(dri 



we have 



dm 



eq(dm-i) 



So, a, e q{dm-i) can be loosely interpreted -the true defini- 
tion being in all cases the relation above- as the epoch of 
equality between the energy-densities of the Dark Matter 
and species i. Its true interpretation is that it provides 
a constant value which measures the ratio of the energy 
density of a relativistic species to the present-day nonrel- 
ativistic Dark Matter energy-density. At times where the 
Dark Matter is relativistic, we can define the ratio 

a eq(dm) = ^ dm^ I ' POdm , 

which we denote just by analogy with a eg ( dm _ i ) as a eq(dm) . 
It is the ratio of the relativistic Dark Matter energy- 
density extrapolated to the present epoch to the actual 
(non-relativistic) energy-density. While, again, the true 
definition is the relation above, when evaluated at a = a nr , 
we see this quantity to be closely related to the offset of the 
energy-density extrapolated from the relativistic side and 
the energy-density extrapolation fron the non-relativistic 
side, that is to the annihilation factor at this epoch. It is 
introduced to provide a measure of this annihilation fac- 
tor, and is of interest only in the NRFO scenario. 

The epoch of decoupling of Dark Matter with species i 
is denoted a dec(dm _i). We have especially a dec(dm ^ v) and 
a dec(dm--y)i the epoch at which Dark Matter decouples 
from neutrinos and photons respectively. The epoch at 
which Dark Matter decouples from all species including 
itself corresponds to a dec ( dm ). We introduce in addition 
the scale-factor which defines the epoch of decoupling of 
neutrinos with electrons, namely a dec ^- e y 

For completeness, we recall that a rec is the scale-factor 
corresponding to recombination. Needless to say, its value 
is taken from the recent WMAP results. The epoch of 
the non-linear collapse, corresponding to a n i, will also be 
needed for strongly interacting Dark Matter particles. 

The analytical and numerical expressions of all of 
these scale-factors are summarized in Table IA.1I The 
Dark Matter properties typically enter these expression 
by means of the Dark Matter mass m dm and interaction 



rate T 



dec(dm-i) 



or r 



dec(dm) • 
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Table A.l. Specific scale-factors . in Sect IA. 2^1 By convention, we define H eq (dm-i) as 

g eg (dw-i) = a eq (d m -i)H r . Together with the rates given 

in Section rA.3.11 this sets the typical scales, displayed in 
anr ~ m dmKdm (T nr ) Table 1X721 to which the Hubble rate H and the reduced 

— 7.06 10~ 10 ( mdm i K & d<r^ r T ™r) ^ interaction rates T of Sect. lA.2~6l are to be compared. 



a. 



eq(-y + v) 



P T + „(T ) 

2.85 10~ 4 



( n m h 7n 2 \ 

V 03 J 



9 ,(T e q)a aq ^ + ^ 



(species i relativistic) 



a eq(dm — i)\ 1 ) — „ 



POdr, 

ig T , z (T) n„ 



O im a eq(f+u ref) 



3.41 10- 4 ^£E1 (^^.y 1 



Pqdm 



-a. 



3.41 10 



3 n 

4 4g'. dm (T) 
3 



eq(-y + f ref) 



( n dm h 70 2 \ 1 

\^ 0.25 J 



Table A. 2. Specific reduced rates . 



H, 



H, 



eq(dm-i) 



NRFO scenario : 

(0dm evaluated at an epoch where it is relativistic) 
(T) 



TtIi 
H„l 



1.48 lO- 29 ^ 1 ,?'? 

in the radiation era 



2.33 10" 32 s 



1 / n m h 70 ' 2 \ 2 ( m dm K dm \-% 

\ 0.3 J \ lMeV I 

atter era 

Vl(T B ,)(a«^Di)" 1 . 

1f1 -24„-i v t< (T) /nwr 1 

iU A 3 » * I 0.25 ^ 

r d ec(,-e) = 2.80 lO- 30 ^ 1 ff'I^-l , 

3 -1 / n m h 70 2 \ 5 

6 i 0.3 y ' 



in the matter era 
5.96 10 _24 s" 

7.16: 



3.45 10" 
cr T hCn e = 5.49 lO"^ 1 s" 



3.93 lO" 20 ^ 1 (ff) 1 (2^)' 



&dec(dm — i) 



O'dec^v — e) 



a r ,_ 



a n I 



2.85 10 







2.85 10 



5.96 10- 24 s-! 

in the radiation era 

?, 



f (T, 



de C (<im)) A. 3. 4. Cross-section scales (Table |A 



(±5-0 
5.96 10- 24 s- : 



2/3 



/ n m h 7 
V 03 



in the matter era 



1.33 10" 



•k 3 



= 1/1090 



1/10 



(but left as a free parameter) 



Our results may also be displayed in terms of the Dark 
Matter (av) cross-sections, more directly connected to 
particle physics. The latter are defined in terms of the 
interaction rate by : 



I\j = {av) i3 rij . 

The values which set the scales of the (av) cross-sections 
relevant to the limits we seek in this work (see e.g. the 
_relations given in Table IB. II which are used in Appendix 



IB.2p turn out to be for neutrinos 



A. 3. 3. Typical Hubble rates fTablelA~2|). 

One can translate the above scale-factors into reduced 
Hubble rates H in order to define typical rates. The re- 
lation between H and the scale-factor a may be found 



S v = ^ - 

while for photons we are led to use the values of the (av) 
-xross-sections related to the epoch of recombination 



or values of (av) related to the Thomson cross-section 
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Finally, we will also be led to use, for any relativistic 
species i, 



H. 



q(dr. 



i) 



4, 

While the definition of s eq is always the above relation, 
it may be (somewhat loosely) interpreted as the (av) 
cross-section needed for the Dark Matter to decouple from 
species i just at the epoch of equality of the energy- 
densities ^ and podm/a 3 - Note that s eq , which may be 
readily evaluated as s eq = Pod Hr /^T uieans of the ex- 
pressions of Heq(dm-i) ( Table IX2)l and rjl i (Section EHl|l ; 
turns out to be independent of the relativistic species i 
considered. 

These typical cross-sections are given in Table IA.3I 



Table A. 3. Specific cross-sections 



Sv =4.48 lCT 1 " 

s rec = 5.53 10~ 26 c 

s Th = 1.99 10" 

s eq = 1.15 10" 



s- 1 9 'i (^my\-i , 

-1 ( rwiToiA -1 

I 0.25 J 



Appendix B: Density evolution and relations 
among interaction rates. 

B.l. Density evolution. 

In case the Dark Matter is non relativistic at the epoch of 
reference, and if the number of particles is conserved till 
nowadays, one has 



ndn 



Qdm POc 



rridn 



In the NRFO scenario, stritly speaking, this relation is 
valid only fr a > af a ~ 7a nr . In the URFO scenario, it is 
valid even for relativistic Dark Matter . 

In the NRFO scenario, for relativistic Dark Matter , 
the number density may be written 11 



ndn 



4 dm K dm 



where n 1 is the density introduced in Section IA.2.51 As 
expected, ndm does not depend on the Dark Matter par- 
ticles' mass. It is independent of time except for possible 

11 The exact relation differs by a numerical coefficient very 
close to unity. We use the above approximation because it joins 
continuously to the non-relativistic expression at a — a nr . This 
turns out to be accurate enough for our purpose and yields sub- 
stantially simpler practical calculations (it avoids in particular 
to distinguish fermionic and bosonic Dark Matter ). 



variations of the factor Kdm(T) or those due to the coef- 
ficient g*dm contained in "4 dm . The densities on the rela- 
tivistic and the non-relativistic side in the NRFO scenario 
are related by 

~ irp rp \ ~i rrp \ K dm(T nr )a nr 

ndm\-L *^ J-nr) — Tdm\ nr > 

^eq(dm) 

Although the Dark Matter density does not depend 
on which species i is present in the system, a quite useful 
relation relating the non-relativistic Dark Matter density 
to the density (see Sect. IA.2~5)l of a relativistic species 
i is 



^ dm 

(T<T n 



4 



(T, 



nr J^nr 



^eq(dm-i) 



B.2. Expression of the I\_d m rate and the {crv) i _ dm 
cross-section in terms of Tdm-i ■ 

The relation (|2T)l between the rates Ti_dm and Tdm-i is 
governed by the ratio of the energy densities i> dm i 1 ^ rela- 
tive to these two species : 



_ f dm r 
-dm — / -L dm— i 



The expressions of Ti_d m may be transformed into expres- 
sions of (^v)i_dm by means of 



/ i — dn 



L i—dm 
T^dm 



To study the correspondence between the interaction rates 
and cross-sections, we need to establish the evolution of 
the number density of Dark Matter particles, that may or 
may not be the same before and after the non-relativistic 
transition. 

We seek a relation between Ti^dm and Tdm-i ■ In the 
special case it is evaluated at an epoch a such as the Dark 
Matter particle number is conserved afterwards, the lat- 
ter can be written in a more explicit form than the gen- 
eral relation 1)21(1. To this purpose, we evaluate the ratio 
fail fadmi specifically under this assumption, for relativistic 
and non-relativistic Dark Matter particles. 

B.2.1. Relativistic species i. 

This is the case we explicitly need to consider in the 
present paper. 

In case the Dark Matter particles are already non- 
relativistic, that is for a > a nr , for both the URFO and 

the NRFO scenario, we have = 



If the Dark Matter particles are still relativistic, on the 
other hand, for a < a nr , we get in the URFO scenario 12 

hm ~ POdmUnr/a 4 , and thus ^ = ■ 



12 The exact relation differs by a numerical coefficient. The 
density ridm in section lFi.ll fsee the footnote there) is related to 
this expression of the energy-density by taking 3Td m as the av- 
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In the NRFO scenario, on the other hand, when a cor- 
respond to an epoch where the Dark Matter still is rela- 
tivistic, since the annihilation is still to come, the Dark 
Matter radiation density may be expected to be of the 
order of the density of any other relativistic species. More 

precisely, for a < a nr , one has — 9 * dm , The ratio 

9 ' *il 1 9 ' *dmi except in very exotic models, is expected (tak- 
ing into account that we are in the NRFO scenario) to be 
of order unity. 

These relation, which will be needed throughout the 
whole paper, and the inferred relation among interaction 
rates and cross-sections are recalled in Table |Bjl 



a quantity adapted to relativistic species i, to use another 
ratio suited to nonrelativistic species i : 



fdm/i 



POdrr 
PiCL 3 



In case the number of particles i is conserved, one has 
obviously f dm/l = 



Table B.2. Relation between i—dm and dm—i interaction 
rates and cross-sections for a nonrelativistic species i. 



Table B.l. Relation between i—dm and dm—i interaction 
rates and cross-sections for a relativistic species i. 



( aV )dm-i = 



a > anr URFO & NRFO 

i>d m _ a 

7 t 



°-eq(dm-i) 



0-iqid.rn.-i 



i — dm an 



a > a nr URFO & NRFO 
- f 

1 Jdm/i 



■Tj— dm — fdm/i ^d 



m d m 1 d m -^ _ m dm I— V 
mi ~ ~ m; \ aU /dm-i 



a < a nr URFO 

^. a e q (dm-i) 



-p Q"nr 

H-dm - a<sq{dm _ t) 



a < a„ r NRFO 



a < a nr URFO 



f) J dm/ 1 a 



^ i — dm 



fdm/i ^ dm — i 



_ Onr m dm r d 
i — dm a 



a < a nr NRFO 



Pdm. _ f , a eq(dm) 
I Jdm/i a 



V. . — g 'dm. p 

-L i — dm — g' • dm — i 



Kdm.{o-v) i 



P _ f a cq(d m ) p 

l i — dm Jdm/i a ^ dm — % 



anr ™-dm 1 dm-i _ a nr . m dm (Tfy) 
a mi a mi \ to 



B.2.2. Nonrelativistic species i. 

Similar expression may be established in this case. They 
will turn out not to be needed in the practical applications. 
We give them here (Table IB. 2(1 for the sake of complete- 
ness: in the remainder of the paper, only the expressions 
for relativistic species i will be given. To this purpose, it 
turns out to be convenient, rather than a e<? (d m _i), which is 

erage energy of a particle independently of the statistics. This 
insures consistency with our choice to set the non-relativistic 
transition at mj m = 3T<i m . 



B.3. The non-relativistic transition in the NRFO 
scenario. 

In the NRFO scenario, there is a sudden decrease of the 
Dark Matter particle number when the latter becomes 
non-relativistic. We approximate the transition around 
a = a nr (and more precisely between a = a nr and 
a = a fo ~ 7a nr ) by an abrupt transition at 
where some physical quantities (such as for instance the 
density) switch abruptly from a "relativistic" value to a 
"non-relativistic" value. This results in a great simplifica- 
tion in displaying our results, and is of no harm, provided 
the associated discontinuities at a = a nr are properly eval- 
uated. 
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It is readily seen from Section lrj.il that in the NRFO 
scenario, there is a drop of the density by a factor 

^dm(^nr) \nr &nr 



Tldm (Tnr ) \ r 



Q>eq(dm) 



This ratio represents the annihilation factor, due to the 
decay of the Dark Matter particles between a nr and a/ OI 
needed to insure cosmological consistency. To represent 
really annihilation it is obviously required to be smaller 
than unity. This makes sense only in Regions I, II, III 
{aeqit+y rcf) > a, nr ) since we expect a eq ( dm ) ~ a eq ( 1+l/ rof) . 

The ratio of the non-relativistic value i> dm (a nr )\ nr to 
the relativistic value /fidm( a nr)\r at a nr is also, from our 
previous estimates (Section LB. 211 of the ratio $J '$ dm 

^dmC^ nr )\ r a eq(dm) 

The question which has been treated in Section IB. 21 is, 
given the interaction rate Td m -i, taken at a fixed epoch, 
what is the expression of the rate Ti-dm and the cross- 
section (ov) i _ dm , taken at the same epoch. We give explic- 
itly here the discontinuities of the latter around a = a nr . 

The rate Td m -i is obviously continuous around a nr . 
From the basic relation H21|) . we then readily see that 

— dm (^rir) | nr Qnr 



-dm 



(T nr ) 1 T 



^eq(dm) 



It bears the same discontinuity than the densities. From 
its definition, (&v)i_ dm = ^i-dm/ndm is the ratio of two 
discontinuous quantities, but remains exactly continuous, 
as it should be, despite the approximations we have made. 
This is possible only thanks to the care we have taken 
in evaluating the approximate expressions of the energy- 
densities, and was by no means granted. 



Appendix C: Classification of the Dark Matter 
candidates. 

As explained in Section [21 the calculation of the damp- 
ing scales is separated into several contributions, the self- 
damping & free-streaming, neutrino and photon induced- 
damping. We treat each of these effects separately in 
order to display the different constraints in a more ap- 
propriate way. The damping lengths turn out to depend 
on the Dark Matter particles' mass and on an appropri- 
ate interaction rate. Their values thus may be displayed 
in a two-parameter space, respectively {rridm, ^dec(dm)]^ 

[mdm,^dec(dm-u)}, [m dm , ^dec(dm-j)] for the self-damping 
& free-streaming, the neutrino induced-damping and pho- 
ton induced-damping. The analytical expression of the 
damping takes different forms in the various "regions" of 
the relevant parameter space, as is seen in Sections 14.11 
15.11 and 16.11 respectively. Here we give the borderlines of 
these regions. 



C.l. Dark Matter classification from self-damping & 
free-streaming. 

The borderlines of the regions defined in Section ^, ll within 
which the damping lengths take different forms are given 
by the conditions 



^dec(dm) 
Q'dec(dm) 



d-nr , 

a eq(7+i/) 
CL nr 



They cross at one single point, and define six regions, la- 
belled from I to VI. 

The discussion in terms of a e9 ( 7+l/ ) is useful because 
this scale-factor remains defined (and relevant) even in 
regions IV, V and VI (although, in this case, it bears no 
relation with the epoch of equality which occurs at a = a nr 
in these regions). 

The numerical values of the above scale-factors are 
given in Appendix I A . 3 . 2| They imply the bordelines given 
explicitely in Table 



Table C.l. Explicit borderlines for self-damping and free- 
streaming 



= 2.48 10~ b MeV ■ 



O'dec(dm) 



& r 



dec(dm) 



5.96 10- 24 s-VJ(T dec(dm) ) 
in the radiation era 



.3 



r. 



dec (dm) 



-i ( n m h 70 2 \ 1 

V 03 J 



in the matter era 



&dec(dm) 

^dec(dm) 



«*• r. 



dec (dm) 



2.33 



in the radiation era 

in-32 -1 ( (t m h 70 2 \ 5 ( m dm K dm y 

iU b \ 0.3 J V lMeV ) 

in the matter era 



C.2. Borderlines of the Regions in the neutrino 
parameter space. 

The damping length due to neutrino induced-damping is 
to be evaluated at Dark Matter-neutrino decoupling. All 
constraints on the interaction rates thus are obtained at 
this time. All relations given in this section hence are un- 
derstood as being written at the epoch a = ct, de c{dm-v) > the 
"reference time", corresponding to the decoupling of the 
Dark Matter with neutrinos. 
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Most important in this case are the interaction rates. Table C.3. Explicit borderlines in the URFO dm — v 
They determine whether v — e oxv — dm scattering or sim- scenario. 

ply neutrino free-streaming is responsible for the damp- 

ing, as discussed in Section [Ql Depending on these rates, _ _ _ 

there are three regions in the [m dm , i'dec(dm—u)] plane, A-B Y dm - V <Y v - d m Y v - d m = Yv- e 

namely A, B and C in which the induced-damping scale 

results from different physical processes : a e q (dm-v) < a nr 

p _/Vgm-»)\3 r 

-L dm — v — I a J ± dec(v-e) 



Region A : r„_ e > r„_ dm , Y dm _ v < Y u ~ T u _ e , Y v - dm = ( a '^ dm -^\ 'r. 



de.c{y — e) 



— Region B : Y u ^ e < T v -& m , Y dm ^ v < ~ T v -& m , 

— Region C : Y dm — v > Y 

Table C.2. Borderlines for dm — v induced-damping. 



A-B borderline 

Ydm — v ^ Yis — dm Yy — dm — Yi/- 

A-C borderline 



&eq(dm — v) CLnr 

( ^dec(dm — iy) ^* a nr 

P _ a d ec (u-e) ^ 

i-v-dfn — a 6 q^d m _ u ^ dec(i/ e) 

/ \ a dec(v-e) 



dm — is -L v — e 



O'dec(dm-u) ^ 0>nr 
I" K u — dm 

Kdm{aV)„_ dm = Su 



p <f P V V v — dm a 1 de.c(y — e) 

J- dm — v _ J- v — dm -L dm — v — -L i^ — e eq(dm. — i/) 

B-C borderline 

r?y— e <c r^—^jTT, rdm-f r^—^jyi 

and (NRFO): r„_ e < r dm _„ < H eq(dm _ v) r dm _„ = H nr B " C r — < r --*n r dm _„ = IV 

a d ec(v — e) ^ ^dec( d m — v) 



dm — v ^-eq( d m — v) 
Yv — dm H eq( d m — v) 



a d ec( d m — v) ^ tt n7* 

The borderlines A-B, A-C, and B-C which separate the p _ 3- 

| _ ^ t t r dm — v ±1 1 

three regions (Table IC.2|) may be directly inferred from 

these definitions. „ "* K *Z~ a ei(dm-u) „ 

Kdm\av) u _ dm — — — — s eq 

It is however of interest to express these relations in 

terms oim dm and Y dec(dm _^. They may also be translated a dec{dm _ v) < a nr 

also in terms of m dm and the cross-sections (&v)„_ dm - We a eq { d m-v) = a„ r 

will do this in the following two subsections (separating 

the URFO and the NRFO scenario) by using the relation 

between Y u _ dm and Y dm _ v established in Appendix IB. 21 D , , .. „ , , s 

vv 1 1 B-C borderline. For a dec ( dm _^ = a nr , that is for Y dm - V — 

H nr , there is a borderline where T v -d m is discontinuous. 

C.2.1. URFO dm-v scenario (Table OJ- On the relativistic side, we are in Region B, with f dm _„ = 

H nr ~ ^is-dm- On the non-relativistic side, we are in 

C.2.2. NRFO dm-v scenario (Table Q. Regi ° n ° With Tdm ^ = Hnr > ^~ dm = °^-* Hnr - 

This however can only be the case for a^ec^-e) < a nr- 

Here, we only consider the cases for which a eq(dm _ u) > Q 3 BordeHines of the Regions in the photon 
a nr . The case corresponding to a eg(rfm- ^) < a nr, is very parameter space 
unnatural as discussed in Section T4. 2. 61 

A-B borderline. For a dec(dm _ !/ ) > a nr , as in the URFO There is now one single borderline since there are only 
scenario, we have a e? ( dm _„) > a nr , which is the case that two regions A and B of interest to compute the photon 
we have decided to not consider explicitely. induced-damping scale. 

All relations given in this section are understood as 
being written at the epoch a = a dec t dm _^\, the "reference 
A-C borderline. For a dec[dm _ u) < a nr , i.e. relativistic U me", corresponding to the decoupling of the Dark Matter 
Dark Matter , one would have ^dm—v — ^u—dm-, which with photons. 
is not possible in Region C. 
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Table C.4. Explicit borderlines for the NRFO dm — v 
scenario (Excluding the case a nr > a eq u m - v \, as discussed 
in the text). 



A-B 



^ dm — v ^ r ^u — dm — dm — Ty- 



^dec{dm — v) ^ Q/nf 



dec(f—e) 



URFO scenario leads to a borderline, given in Table 
where we have Tdm-~f = eg( a dT "~ T) r Th . For all sensible val- 
ues of a nr namely a nr < a eq ^ d 

m— f^Th/ Heq(dm— 7) ! it corre- 
sponds also to Dark Matter which decouplesjn the matter 
dominated era. It is also readily seen that F<jm— 7 < H rec 

implies a nr > a e q(dm A^t^I Hrec- Note that the limit is 

reached at a value a nr = a eg (d m - 7 )r Th /iJ ree > a e q (dm-~/)- 
The second borderline (see Tables and RXYjl . in the 
URFO as well as the NRFO case, corresponds to 



A-C 



and is relevant for the smaller values of a nri namely a nr < 
a eq(dm--y)^ Th / Hrec, that is for the larger masses rridm- 



^dec^dm — v) ^* a vir 
^dm — u ■ ^dec(v-e) 



IV 



a dec(v-e) 
a eq(dm-v) 



^ dec(u — e) 



/ \ "dec! 

Kdm{arv) v _ dm = —± 



Table C.6. Explicit borderlines for the URFO dm — 7 



scenario. 



A-B 



f^dm — 7 ^ Hrec 



- 7 — dm 



B-C 



^dec(v-e) ^ ^ , dec(dm — u) 

/ \ a eg(dm-y) 



&dec(dm — 7) a rec 
■ dm — 7 — 
1^7 — dm TtIi 



P _ a eq(dm.- T ) ^ 
J- dm — -/ — ~ 1 Th 



Kdm{(Tv} y _ dm = K — 



B-C 



B&C 



eq(dm-u) -L dm—u 



Q>dec{y — e) ^■dec(dm-u) — Q>nr < a eq(dm — v) 



B side: (where T v - e < IV ~ T 



v — dm ~ J- dm — v ) 



IV 



C side: (where r«_ e < r„ ~ IV dm < r dm _„) 



r„ 



1 Hnr 



a eq(dn 



A-B 



( ^ J dec(dm — 7) ^* 



( ^dec(dm — 7) ^ fl nr 



- dm — 7 ■ — H rec 



^dec{dm — "{) — &rec 



p a dec(dm- T ) 77 

J- 7 — dm „ , , „ "f> 



Kdm{<JV) 



a dec(dTn.-7) 



7 — dm 



7 — dm — „ J^rec 
S(dm-l) 



7 — dm 



Table C.5. Borderlines for the dm - 7 scenario. Table C.7. Explicit borderlines for the NRFO dm - 7 

scenario (Evaluated only for a„ r < a eq ^ m ^j as discussed 
in the text). 



A-B 


Tdm — 7 *C H rec 


1^7 — dm 


— r 

J. -y — e 


A-B 




Tdm — 7 





A-B r,j m _ 7 = _ff 



This may be rewritten in terms of r,j m -7 by means of the 
relations of Section IB. 21 Relativistic Dark Matter , in the 



^dec(dm — 7) — &rec 



Provided Tdrn-~t < H rec , the first borderline (Table 
IC.5(1 corresponds to p _ a d SC (d m -~i) jj 



J- 'V— R 1 Tb ■ ' 



"d e c(dm- 7 ) 



7 — dm a 711 - ^rec 
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Appendix D: Damping limits. 

We give in this section the analytical form of the constraint 
one obtains by writing that the various damping lengths 
have to be smaller than a given a priori chosen limit I struct ■ 
Note that in each case, the parameter space is constrained 
by the relevant damping limit and the condition to be in 
a given "Region" in which the expression used for the 
damping scale is valid. 



D.l. Self-damping. 

The damping scales in this case are given as a fraction 

i 

of the horizon at decoupling 7rct r a dec(dm) , or ^ct m a^ ec(dm) . 

i 

Noting that t m differs from t r basically by a factor ai q 
(Sect. IA.2.1|) . one sees all damping lengths can be ex- 
pressed as a fraction of 



irct r a e 



207 Mpc g' 



0.3 



The analytic expression of the damping scales and the as- 
sociated limits on the interaction rates and the (av) cross- 
sections are given in table ID. II for the six Regions which 
cover the parameter space. The free-streaming damping 
scale has been calculated via the approximate expression 
<|»T7|) . to be directly comparable to the expression of the 
self-damping length, estimated using a similar approxi- 
mation (|29[) and to be related to the mass scale by l|38|) . 
In Regions I, II and III where our limits are set, we give 
for completeness -within brackets- also the expression of 
the free-streaming scale when it is by convention taken to 
be given by (|36|l . All relations given in this section are 
understood as being written at the epoch a = ad ec (dm), the 
"reference time", corresponding to the decoupling of Dark 
Matter with every species including itself. 



Table D.l. Self-damping limits 



Region I 

lfs 

[ hs 
Isd 



= L 



L eqr dm 

< a eq — 



5+21n(a cq /a Ilr .) 



Region II 

l f s ~ r rm ~ e< ?l 
-/ (- 



Ifs 



Region III 

[ hs 



Region IIP 



r a dec(dn 



1 3 + 21n(a e q/q dee ) 



Region IV 

lfs 
hd 



Region V 

lfs 
Isd 



a dec(d; 



— Lteq I , 



1/2 



= L 



'eq^ dm 

< a eq — 



<■{<<•••) 



= I 
= L 



<:q , 



l &qr d 

I struct 



a Jec(Jm) 



D.2. Neutrino induced-damping. 

The typical scale in this case is the damping scale of the 
neutrinos at the time they decouple from the electrons 



Region VI 
l fs = 7± 



<aeq (k^y 



ldec(v—e) — 7TCt r ri/- 



dec(v-e) 



97.1 pc Tv K 



For Region A, the damping is too small to be rele- 
vant to structure formation. For Region B, only the case 
a dec(dm~v) < a eq(dm-v) has been considered explicitly: for 
adec{dm-u) > a e q{dm-v)-, the damping is too large. Needless 
to say, in the latter case, too, the damping may be readily 
evaluated by means of the same methods. 

Analytical expression of the damping scales and limits 
on interaction rates and (av) are given in tables lD~2l ID. 31 
ID. 41 ID. 51 and, ID. 61 All relations given in this section are 
understood as being written at the epoch a = a>dec(dm-v)> 
the "reference time", corresponding to the decoupling of 
the Dark Matter with photons. 



Region VI' 

!.,rf _ / ( 3 (a nr \l/2 



D.3. Photon induced-damping. 

A unit to characterize the damping scale in this case may 
be chosen as the horizon scale at equality, whcighted by 
a slightly time-dependent factor r 7 , to be taken at the 
proper time: 



leq(dm-y) 7T^r^7^eg(rfm— 7) 



249 Mpc r 7 g 



0.25 
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Table D.6. Region C (NRFO dm - v scenario). 
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^dec(u-e) {J^drn — h'/^dec(v — e) ) * 



Table D.3. Region B (URFO dm - v scenario) 



— tdec(^-e) I J— )■ 



1 dcc(v-c) 



< r 



htr 



Table D.4. Region B (NRFO dm — v scenario) 



< r 



7 ( J™ \ * T *™-» 

Struct 



(»-*) 



dec(y — e) ' 



- dm-v 



< r 



dec(u-e) i 



I struct 



dec{u-e) 



&dec(drn — v) ^* Q>nr 



a dec(dm,-is) T drl 



Q>dec(dm — iy) ^ a nr 

does not exist in this case 



matter dominated era (a,dec{dm-~() > a eg)- We do not give 
the damping lengths when Dark Matter decouples from 
the photons in the matter dominated era since they are 
always prohibitive as compared to our requirements. 

Analytical expression of the damping scales and limits 
on interaction rates and (av) are given in Table ID. 71 for 
decoupling in the radiation dominated era. They are ob- 
tained from the expression of the damping length given in 
Section refsec:photlim by means of the relations given in 
Table IB~T1 and the definitions of Appendix lAl All relations 
given in this section are understood as being written at the 
epoch a = 0^(^-7). the "reference time", corresponding 
to the decoupling of the Dark Matter with photons. 



Table D.5. Region C (URFO dm — v scenario). 



Table D.7. Region A, radiation dominated (URFO and 
NRFO dm — 7 scenario). 
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< S 



°fa(i-0 f htruct V 



Kdm(crw)„_ 



< Si 



dm — j 



eq(dm — 'y) 



&dec(dm — 7) ^ Q nr 

a dec((im-7) T dm _ 7 
— — r J~ 



K dm {av) dm = 



< S e 



( htruct \ ' 
\ l cq{d m -^) ) 



eq(dm--,) ( htruct \ 
\ l cq(dm--,) ) 



In Region A, rather than the URFO and NRFO sce- 
narios which turn out here to yield the same damping 
limits, it turns out to be useful to distinguish two cases, 
depending on whether Dark Matter decoupling occurs in 
the radiation dominated era (<idec(dm~^) < a eq), or in the 



a, 1,. 



■(dm — 7) 
■{<™)j-dm 



r rf ,„ 



< S e 



I htruct 
\ l cq(d m —, 



